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1.   Introduction 

This  document  presents  results  from  the  Lompoc  Valley 
Diffusion  Experiment  (LVDE) ,  completing  tasks  related  to  the  LVDE 
data  analysis  described  in  contract  number  MPIRF7616890425  funded 
by  USAF  HQ  Space  Systems  Division,  Los  Angeles  (SSD/CLGR) .  This 
report  uses  data  sets  described  in  the  LVDE  Data  Report 
(Skupniewicz  et  al.  1991a),  and  therefore  detailed  descriptions  of 
the  measurements  and  methods  of  LVDE  are  not  supplied  herein.  The 
LVDE  Data  Report  also  supplies  some  basic  analyses  of  the  data, 
such  as  plume  trajectories,  regional  wind  flow  patterns,  and 
synoptic  descriptions. 


The  first  objective  of  this  analysis  is  to  compare  LVDE 
results  to  the  Mt.  Iron  equations,  and  to  formulate  new  expressions 
specific  to  LVDE.  During  this  procedure,  we  simulate  the  methods 
used  in  the  Mt.  Iron  reports  (Hinds  and  Nickola  1968)  .  We  discuss 
differences  in  the  results,  and  show,  some  difficulties  with  the 
statistical  approaches  used  in  the  Mt.  Iron  analysis.  The  second 
objective  is  to  compare  LVDE  data  to  theories  developed  for 
diffusion  of  passive  materials  from  surface  sources  in  convective 
conditions.  We  extend  those  basic  theories  to  accommodate  the 
added  complexities  of  the  LVDE  domain. 


2 .   Background 

For  those  readers  without  access  to  the  LVDE  Data  or  Mt .  Iron 
reports,  we  describe  these  experiments  briefly  in  this  section. 
LVDE  was  conducted  during  August  1989  to  assess  the  potential  risk 
of  gaseous  releases  from  the  Hypergolic  Stockpile  and  Storage 
Facility  (HSSF)  at  Vandenberg  AFB,  California.  Previous  studies  of 
wind  flow  patterns  and  hypothetical  gaseous  releases  showed  that 
trajectories  could  pass  over  populated  areas  near  Lompoc, 
California  during  typical  daytime  onshore  wind  flow  (Kamada  et  al. 
1989)  .  Thus,  tests  were  conducted  under  normal  daytime  "sea  breeze" 
conditions.  Releases  were  near-surface  (12  ft)  and  continuous. 
The  trace  gas  was  SFs  with  downwind  sampling  by  mobile  electron- 
capture  gaseous  samplers.    The  instrumented  vehicles  travelled 


roads  which  intersected  the  plume  at  2-15  km  ranges.  LVDE  also 
included  a  series  of  releases  from  within  the  adjacent  Lompoc 
Valley  2  km  north  of  the  HSSF.  We  term  these  "flat  terrain"  data 
as  opposed  to  the  hilly  complex  terrain  downwind  of  the  HSSF. 

The  Mountain  Iron  Diffusion  Tests  were  a  series  of  surface 
releases  conducted  in  the  1960's  from  Space  Launch  Complexes  4  and 
6  (SLC4,6).  The  original  report  developed  a  range  of  regression- 
based  equations,  depending  upon  the  release  site  of  interest  and 
the  input  measurements  available.  This  report  uses  a  concentration 
vs.  range  equation  developed  for  SLC4  which  is  2  miles  south  of  the 
HSSF.  Terrain  downwind  of  SLC4  is  similar  in  complexity  to  the 
HSSF  downwind  terrain.  The  releases  were  typically  30  minutes  in 
duration.  The  trace  material  was  a  fluorescent  aerosol  (zinc 
sulfide,  GMD  2.5  microns) .  Downwind  sampling  was  accomplished  with 
vacuum  samplers  distributed  along  available  roads  on  the  base.  The 
results  from  Mt.  Iron  are  presently  used  for  assessment  of  neutral 
gas  surface  spills  from  locations  in  this  general  area. 

While  LVDE  and  Mt .  Iron  were  conducted  under  similar  wind  flow 
patterns  and  the  release  locations  were  within  4  km  of  each  other, 
some  fundamental  differences  in  the  experimental  designs  are  noted 
here.  First,  Mt .  Iron  released  tracer  for  relatively  short  periods 
of  time,  typically  30  minutes,  while  the  LVDE  release  durations  of 
four  to  seven  hours  were  essentially  continuous  (when  compared  to 
the  travel  time  of  the  plume  to  sample  locations).   Secondly,  the 


Mt .  Iron  trace  material  was  collected  over  the  duration  of  the 
release,  giving  a  dosage,  or  time  integrated  concentration,  while 
LVDE  data  were  gathered  from  rapid  transects  through  the  plume 
giving  "instantaneous"  crosswind  concentration  profiles.  Plume 
parameters  calculated  from  the  LVDE  data  can  be  modeled  with  two- 
particle  statistical  theory,  useful  in  estimating  the  crosswind 
growth  of  instantaneous  releases  ("puffs").  When  time-averaged, 
the  LVDE  data  simulates  the  growth  of  continuous  plumes,  described 
by  single-particle  diffusion  theories.  The  time-averaged  form  of 
the  LVDE  data  is  more  comparable  to  the  Mt.  Iron  results  than  the 
instantaneous  data. 


3.   Mt  Iron  Regression  Analysis  Methodology 

Risk  assessments  of  most  types  of  gaseous  releases  at 
Vandenberg  are  primarily  concerned  with  the  hourly  averaged 
surface-level  maximum  concentration.  The  concentration  is 
integrated  over  the  time  of  exposure,  resulting  in  accumulated 
dosage.  The  USAF  Surgeon  General  has  set  maximum  dosage  levels  for 
various  gases  based  on  toxicology  studies  and  standards  set  by 
other  public  health  organizations  (e.g.  OSHA) . 

A  primary  goal  of  Mt.  Iron  was  to  predict  centerline 
concentration  using  measurements  from  meteorological  towers  located 
at  or  near  the  release  sites.   The  tower  data  is  readily  available 


at  Vandenberg  in  near  real  time.    Investigators  chose  a  multi- 
parameter regressive  approach,  assuming  that 


-5?  =  f   (x,oe,AT,U)    =  a^ap^T^U^     ,  (1) 


where  Cm  is  maximum  ground-level  concentration,  Q  is  release  rate, 
x  is  downwind  range,  oq  is  standard  deviation  of  the  wind 
direction,  AT  is  temperature  difference  between  tower  sensors  at  54 
ft  and  6  ft,  and  U  is  wind  speed.  The  left  hand  side  of  the 
relationship  is  dimensioned  (ppm  N02~min)/lb.  Concentrations  of  any 
other  chemical  species  are  obtained  by  multiplying  the  expression 
by  the  ratio  of  their  molecular  weights  to  that  of  NO2 .  From  the 
right  hand  side  uses  the  non-conventional  units:  ft,  degrees,  deg 
F,  and  knots.  While  these  units  are  awkward  by  present  standards, 
we  retain  this  convention  to  compare  LVDE  with  Mt.  Iron  directly. 
The  Mt.  Iron  experiment  did  not  measure  concentrations,  but  rather 
dosage  accumulated  over  the  period  of  release.  Results  were  later 
converted  to  concentration  by  assuming  uniform  concentration  over 
the  release  period,  and  dividing  by  the  release  period  and  sampler 
flow  rate. 

The  regression  parameters  used  in  Mt .  Iron  were  chosen  for 
their  significance  to  the  physics  of  diffusion.    Turbulence 
intensity,   i.e.   cjq  is  of  primary  importance.    The  spectral 


characteristics  of  turbulence  determine  changes  in  diffusion  as  a 
function  of  plume  travel  time,  and  this  can  often  be  related  to 
atmospheric  stability,  e.g.  the  bulk  Richardson  number, 


Rb=MZl 


T 


(    dT 
oz 

\      u2     ) 


(2) 


therefore,  AT  and  U  where  chosen  as  regression  parameters.  The 
centerline  concentration  is  also  a  function  of  the  plume  material, 
release   rate,   distance   from  the   source,   and  possibly   other 
meteorological  or  physical  quantities. 

It  is  apparent  that  several  of  the  regression  parameters  are 
interdependent.  This  implies  that  multivariate  regression 
equations  formed  by  combinations  of  these  variables  may  not  be 
unique.  Therefore,  it  makes  more  physical  sense  to  group  the 
parameters  in  non-dimensional  forms,  thus  reducing  the  number  of 
empirically  derived  constants.  We  also  note  that  the  Mt.  Iron 
regression  neglects  some  important  physical  quantities,  e.g.,  the 
mixing  depth,  h.  We  will  discuss  these  issues  further  in  following 
the  sections. 


The  original  technique  for  evaluating  the  Mt.  Iron  equations 
was  to  count  the  number  of  occurrences  of  samples  which  fell  within 
factors  of  two  and  four  of  the  equation  predictions.  While  other 
grading  schemes  may  be  more  insightful,  we  have  kept  this  "factor 
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grade"  approach  for  comparison.  We  have  also  evaluated  the  results 
by  using  the  normalized  mean-squared-error  (NMSE)  and  fractional 
bias  (FB)  ,  measures  favored  by  other  investigators  (e.g.  Hanna  and 
Strimatis  (1991).   These  are  defined  as 


(C  -C  )  2 
NMSE  =  p      °  ,  (3) 

C  C 


and 


c  -c 

FB   =  2  _Z   JL      ,  (4) 


where  subscript  o  refers  to  the  observations  and  p  refers  to  the 
equation  predictions.  By  dividing  by  the  mean  observed  and 
predicted  values,  NMSE  tends  to  give  equal  weight  to  under- 
predictions  and  over-predictions.  For  example,  in  the  case  of  two 
observations  of  say,  5.0,  then  a  prediction  of  2.5  has  the  same 
NMSE  as  a  prediction  of  10.0.  In  a  similar  manner,  FB  gives  the 
percent  of  over  or  under-prediction.  Both  measures  suffer  from  the 
problem  that  the  error  characteristics  cannot  be  a  function  of  the 
absolute  value  of  the  data.  I.e.,  the  NMSE  or  FB  at  1  ppm  should 
be  the  same  as  that  at  1  ppb.  If  not,  the  data  at  1  ppm  will  be 
given  more  weight.  We  found  concentration  independence  to  be  true 
for  our  application  of  NMSE,  but  untrue  for  FB  applications. 
Therefore  we  define  an  alternative  measure,  "geometric"  fractional 


bias, 


log  (CJCJ 
GFB  =   °     p   —  ,  (5) 

abs   (log  (C0CD)  ) 


'o~p' 


which  is  similar  to  FB  but  with  a  logarithmic  weighting.  GFB  at 
low  concentrations  are  equivalenced  with  bias  at  higher  absolute 
concentration.  The  disadvantage  of  GFB  is  that  it  cannot  be 
directly  related  to  a  fraction  of  over  or  under-prediction. 

Mt.  Iron  results  are  given  both  in  terms  of  median  value  and 
"95%  confidence"  equations.  We  are  unclear  as  to  how  the  original 
investigators  derived  the  95%  confidence  equations,  but  the  two 
forms  are  related  by  a  simple  fraction,  i.e.,  only  parameter  a^  in 
eq.  1  differs  between  the  two  equations.  We  obtain  a  95% 
confidence  equation  by  assuming  a  normal  distribution  with  a 
variance  of  NMSE.  Then  approximating  95%  confidence  at  "two 
sigma" , 


-^(95%)  =  (1  +  2  /NMSE)     -^(50%)   .  (6) 


Unfortunately,  base  communications  to  the  meteorological  tower 


10 


at  the  release  site  (WT057)  failed  just  prior  to  the  start  of  LVDE. 
While  we  were  able  to  capture  some  data  by  interfacing  to  the 
outgoing  1  minute  data  stream,  these  data  are  not  processed  through 
the  base  computer  system  into  the  standard  10  minute  averages  and 
have  numerous  data  gaps.  Rather  than  deciphering  and  simulating 
the  processing,  and  deleting  data  during  the  met  data  lapses,  we 
decided  to  use  a  nearby  tower  (WT054)  for  the  met  inputs.  The 
tower  chosen  was  about  1  mile  from  the  release  site,  roughly  the 
same  distance  from  the  coast,  in  similar  terrain  with  respect  to 
surface  roughness,  and  under  similar  cloud  cover  as  during  the 
releases.  Met  data  from  two  other  towers  were  also  tried  in  the 
analysis  (WT058  and  WT300) . 

In  Skupniewicz  et  al.  (1991a),  we  calculated  the  centerline 
maximum  concentration  via  three  methods.  The  "moment"  method 
performs  moment  statistics  on  the  cross-wind  concentration 
distributions.  The  other  two  methods  assume  a  Gaussian  crosswind 
profile.  The  "maximum"  method  assumes  that  the  maximum  observed 
concentration  is  true,  then  forces  ay  to  agree  with  the  cross-wind 
integrated  concentration  (CWIC) .  The  last  method  is  based  solely 
on  CWIC,  forcing  both  the  maximum  and  ay  to  conform  to  the  shape  of 
the  profile.  We  use  only  the  "maximum"  method  for  these 
regressions,  since  we  want  the  most  direct  measure  of  maximum 
centerline  concentration. 

Hourly  averages  were  calculated  with  the  following  equations; 
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where  subscript  i  denotes  instantaneous  transects  and  yj_  is  the 
cross-wind  position  of  the  instantaneous  maximum  concentration. 


4.   Mt.  Iron  Regression  Results 

Figures  1-4  compare  the  range  dependence  of  the  LVDE  data  to 
that  of  Mt.  Iron  for  both  instantaneous  and  time  averaged,  and  both 
flat  and  complex  terrain  data  sets.  The  following  generalities 
apply  to  all  figures: 

*   Any  data  points  overlaid  with  "star"  symbols  were  edited 
prior  to  the  regression  analysis,  but  were  included  in  the  grading 
schemes. 
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*  Unless  otherwise  indicated,  the  Mt.  Iron  equation  is  the 
95%  confidence  version  (operationally  used  at  Vandenberg) ,  while 
the  LVDE  data  are  median  (50%)  values.  The  Mt .  Iron  95%  line  is 
approximately  a  factor  of  3  higher  than  the  Mt .  Iron  median  line. 
We  directly  compare  95%  values  later  in  the  report. 
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Figure  1.   Instantaneous  maximum  concentrations  normalized  by  the 
release  rate  vs.  downwind  distance  for  the  HSSF  release  site. 
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Figure  2.   One  hour  averaged  maximum  concentrations  normalized  by 
the  release  rate  vs.  downwind  distance  for  the  HSSF  release  site. 
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Figure  3.  Instantaneous  maximum  concentrations  normalized  by  the 
release  rate  vs.  downwind  distance  for  the  Lompoc  Valley  release 
site. 


C  =  moxi    tower  =  054 


-  ao'tea    line    =     Mt.    Iron    fil 
.     sono    line    =   .'JVOZ    fit 

-  circies    =    LVDC    actio 


C/Q=0.135614X 


-0.591 


: 


10' 


X    (ft) 


- 


Figure  4.  One  hour  averaged  maximum  concentrations  normalized  by 
the  release  rate  vs.  downwind  distance  for  the  Lompoc  Valley 
release  site. 


For  ranges  where  plume  dimensions  are  small  compared  to  the 
dominant  turbulence  scale,  theory  predicts  that  concentrations 
resulting  from  a  puff  (instantaneous  data)  would  drop  less  rapidly 
with  range  than  plume  concentrations  (time  averaged  data) .  We  see 
this  effect  by  comparing  the  instantaneous  data  of  figs.  1-3  to  the 
time  averaged  results  of  figs.  2-4.  Eguations  obtained  from 
curvilinear  range  dependent  fits  to  the  LVDE  data  are  at  the  top  of 
each  figure.  The  exponents  in  the  regression  equations  indicate 
the  rate  of  concentration  decrease  with  range.  Comparing  the 
instantaneous  to  time  averaged  exponents,  we  see  smaller 
exponential  decreases  after  averaging  is  applied. 

In  all  cases,  the  best  fit  through  the  LVDE  data  is  less 
sloped  than  the  line  obtained  by  regressing  the  Mt.  Iron  equation 
against  the  downwind  distance.  This  shows  that  the  Mt.  Iron 
concentrations  decrease  much  more  rapidly  with  range  than  the  LVDE 
data.  This  feature  cannot  be  attributed  to  the  "short"  Mt.  Iron 
release  time  (30  min) .  The  short  releases  of  Mt.  Iron  should  cause 
the  data  to  behave  in  a  more  "plume-like"  fashion  at  short  ranges 
and  more  "puff-like"  at  long  range,  resulting  in  a  greater  slope 
than  a  true  continuous  release  would  show.  This  is  not  observed. 
Several  factors  could  produce  these  differences.  A  fall-out  of  Mt. 
Iron  tracer  material  with  increasing  range  would  produce  this 
effect.  Enhanced  diffusion  at  short  range  (compared  to  Mt.  Iron) , 
or  decreased  diffusion  at  long  range  could  also  produce  these 
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results.  Unfortunately,  regression  techniques  such  as  these  do  not 
allow  for  a  unique  diagnosis  of  these  differences. 

Figures  5-8  apply  a  multi-variate  curvilinear  regression  to 
the  full  set  of  parameters  listed  in  eq.  1,  and  plot  that  solution 
against  the  measurements.  The  data  scatter  within  the  factor  two 
and  four  lines  is  very  similar  to  the  Mt.  Iron  results;  73%  of  the 
Mt.  Iron  data  were  within  a  factor  of  two  of  the  predictions  and 
97%  were  within  a  factor  4.  The  NMSE  values  are  quite  small 
compared  to  other  model  evaluations  such  as  Hanna  and  Strimatis 
(1991),  but  this  is  misleading  because  our  regression  results  are 
being  tested  against  the  data  from  which  they  were  obtained.  The 
improvements  in  the  factor  grades  and  NMSE  after  averaging  are  also 
somewhat  artificial  because  the  number  of  data  points  has  been 
reduced,  and  the  data  sets  are  not  independent.  If  we  employ  the 
central  limit  theorem,  assuming  samples  from  the  same  normally 
distributed  population,  we  would  expect  roughly  a  factor  of  3.7 
reduction  in  the  variance  for  both  the  HSSF  and  Lompoc  Valley 
averaging  processes.  If  we  take  the  ratio  of  the  NMSE  values,  we 
find  reductions  of  2.4  for  HSSF  releases  and  2.1  for  the  Lompoc 
Valley  releases,  well  below  the  central  limit  for  both  cases.  This 
shows  that  the  averaging  process  has  actually  added  some  variance 
to  the  data. 
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Figure  5.  Scatter  plot  of  instantaneous  observed  concentrations 
vs.  multi-parameter  regression  equation  (listed  at  top)  for  HSSF 
releases.  Diagonal  lines  represent  values  within  a  factor  of  two 
and  four  of  agreement.  No  data  have  been  edited  when  calculating 
these  "goodness  of  fit"  grades. 
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Figure    6.        Same    as    figure    5,     except    data    are    one    hour    averaged 
concentrations    from  HSSF   releases. 


-0.932        -1  078  1.185     0.02} 

C/0  =  6.149X  <7  (AT+9)  U 


h  

r     65    %  within    toclor    ol    2 

r    85    r.  within    toclor    of    4 

norm,  mean    sq.    error    0.80 


/  o 


T3 
0) 


O 


3      4    5   5  7  39 


10' 


10 


-1 


C/Q    fitted 


C  =  moxx    io~»r  =  054 


5  7  39 


10' 


Figure  7.    Same  as  figure  5,   except  data  are  instantaneous 
concentrations  from  Lompoc  Valley  releases. 
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Figure  8.    Same  as  figure  5,  except  data  are  time  averaged 
concentrations  from  Lompoc  Valley  releases. 


We  next  discuss  concentration  dependence  on  oq,  AT,  and  U  as 
measured  at  the  local  met  tower.  The  most  disturbing  feature  of 
the  oq  regressions  were  weak  C/Q  correlations  with  measured 
turbulence  for  the  HSSF  release  site  (e.g.,  fig.  9).  We  would 
expect  a  well  defined  negative  relationship  between  C/Q  and  sigma 
theta,  since  higher  turbulence  disperses  the  cloud  more  rapidly. 
The  negative  correlation  is  better  for  Lompoc  valley  releases  (fig. 
10) ,  but  the  scatter  is  large.  The  poor  correlation  is  contrary  to 
most  theoretical  descriptions  of  dispersion,  and  probably 
indicates  a  flaw  in  the  measurement  system.  We  have  determined 
that  the  base  computer  system  does  simple  averaging  of  one  minute 
oq  values  to  determine  longer  averages.  This  high  pass  filtering 
incorrectly  eliminates  much  of  the  turbulence  responsible  for 
diffusing  the  plume.  Steps  are  under  way  to  correct  this  situation 
at  the  base. 

When  we  applied  a  multi-variate  curvilinear  fit  to  the  range 
and  sigma  theta  parameters,  we  obtained  only  a  minimal  improvement 
in  the  factor  grades  and  NMSE  over  the  range  fits  alone.  This 
small  degree  of  increased  predictability  is  largely  due  to  the  lack 
of  C/Q  correlation  with  sigma  theta.  Therefore,  the  exponent  for 
oq  in  figs.  5-8  are  poor  estimates,  and  can  suggest  a  stronger 
relationship  than  truly  exists  (e.g.  fig.  7) .  Indeed,  the  Mt.  Iron 
exponent  may  be  just  as  unreliable  as  the  LVDE  exponents,  but 
without  the  original  data,  we  cannot  confirm  this  hypothesis. 
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Figure  9.   Instantaneous  maximum  concentrations  normalized  by  the 
release  rate  vs.  turbulence  (aQ)  for  the  HSSF  release  site. 
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Figure  10.   Instantaneous  maximum  concentrations  normalized  by  the 
release  rate  vs.  turbulence  (aa)  for  the  Lompoc  Valley  releases. 
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Figure  11.  Instantaneous  maximum  concentrations  normalized  by  the 
release  rate  vs.  vertical  temperature  difference  at  tower  054  for 
the  HSSF  release  site. 
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Figure  12.  Instantaneous  maximum  concentrations  normalized  by  the 
release  rate  vs.  vertical  temperature  difference  at  tower  058  for 
the  HSSF  release  site. 


The  temperature  difference  between  54  ft  and  6  ft  is  a  rough 
measure  of  atmospheric  stability;  more  negative  values  indicate 
more  convection  and  therefore  greater  dispersion.  As  in  the 
original  Mt.  Iron  reports,  our  analyses  add  9  deg  F  to  AT  in  order 
to  stabilize  the  regression  by  keeping  values  positive.  We  expect 
a  positive  relationship  with  AT.  We  found  this  to  be  true,  but 
again,  the  scatter  is  large.  The  dependence  is  weak  for  WT054 
data  (fig.  11)  because  it  was  mainly  under  the  clouds  where 
stability  is  near  neutral.  The  same  data  plotted  against  WT058 
temperature  differences  (fig.  12)  showed  significantly  better 
correlation,  especially  for  low  (more  negative)  temperature 
differences.  AT  is  generally  smaller  at  WT058  because  it  was 
mainly  in  clear  skies  and  unstable  conditions.  Performing  the 
multiple  variable  regression  reflects  this  dependence,  with 
significant  error  reductions  for  WT058  met  data  compared  to  WT054 
data.  A  subtle  conclusion  may  be  made  from  this  exercise.  In  some 
cases,  the  release  site  tower  may  not  be  the  best  choice  for 
calculating  dispersion.  For  releases  under  cloud  cover,  the 
release  site  would  most  likely  be  in  nearly  neutral  conditions. 
However,  the  dispersion  along  the  plume  pathway  takes  place 
primarily  in  clear  skies.  Thus,  better  correlation  may  observed 
for  AT's  measured  in  the  clear  skies,  as  was  the  case  with  WT054 
and  WT058. 

Wind  speed  is  a  natural  dispersion  parameter  for  a  continuous 

29 


plume  because  the  amount  of  material  released  into  a  unit  volume  of 
air  (for  a  constant  mass  release  rate)  decreases  linearly  with  wind 
speed.  When  we  regress  concentration  against  wind  speed  alone,  we 
find  a  strong  negative  relationship  for  all  towers  (figs.  13-14). 
However,  when  we  included  wind  speed  in  a  multi-variate  regression 
with  the  other  parameters,  we  did  not  observe  a  significant  error 
reduction.  This  is  reflected  in  the  wide  variation  of  exponents 
given  in  figs.  5-8.  We  attribute  this  parodox  to  interdependencies 
between  regression  parameters.  For  example,  the  lower  effective 
release  rate  with  increasing  wind  speed  is  compensated  to  some 
degree  by  lower  turbulence  intensity  values,  and  therefore  less 
dispersion. 
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Figure  13.   One  hour  averaged  maximum  concentrations  normalized  by 
the  release  rate  vs.  wind  speed  for  the  HSSF  releases. 
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Figure  14.   One  hour  averaged  maximum  concentrations  normalized  by 
the  release  rate  vs.  wind  speed  for  the  Lompoc  Valley  releases. 


Figures  15-18  apply  eqn.  6  to  the  LVDE  data,  giving  95% 
confidence  figures.  The  recommended  95%  confidence  equation  for 
time  averaged  HSSF  releases  is  given  atop  fig.  16.  These  figures 
show,  in  an  absolute  sense,  the  difference  between  the 
measurements.  For  all  cases  we  see  a  bias  change  with 
concentration.  Figure  16,  the  case  which  should  agree  most 
closely,  shows  significantly  larger  LVDE  values  than  Mt.  Iron 
concentrations  at  low  absolute  values  (longer  ranges).  The  bias  is 
reduced  or  eliminated  at  higher  concentrations  (shorter  ranges) . 
Again,  this  changing  bias  could  be  explained  by  a  fall-out  of  Mt . 
Iron  tracer  material  with  increasing  range,  enhanced  diffusion  at 
short  range  (compared  to  Mt .  Iron) ,  or  decreased  diffusion  at  long 
range.  Unfortunately,  we  cannot  deduce  the  cause  from  these 
regression  analyses. 

Figure  16  demonstrates  that  GFB  may  be  more  appropriate  than 
FB  in  those  cases  in  which  the  bias  changes  with  concentration. 
The  large  negative  fractional  bias  occurs  because  a  few  high-valued 
"predictions"  (Mt.  Iron)  are  greater  than  the  "observations"  (C/Q 
measured) ,  when  in  reality  the  observations  are  generally  larger. 
The  geometric  fractional  bias  is  close  to  negligible,  since  the 
higher  LVDE  values  for  lower  concentrations  are  given  equal  weight. 

Inspection  of  figs.  16  and  18  also  indicates  that  the  time 
averaged  "flat  terrain"  concentrations  are  generally  lower  than  the 
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"complex  terrain"  data,  while  figs.  15  and  17  show  negligible 
difference  for  the  instantaneous  data.  Comparing  the  HSSF  and 
Lompoc  Valley  releases  for  distance  only  regression  (figs.  1-4), 
we  see  that  the  flat  and  complex  terrain  data  are  similar  in  an 
instantaneous  sense,  but  the  flat  terrain  concentrations  (Lompoc 
Valley)  decrease  less  rapidly  in  a  time  averaged  sense.  The  large 
scatter  and  small  number  of  data  points  leave  us  wary  of  suggesting 
physical  processes  for  these  differences.  However,  we  further 
analyze  Lompoc  Valley  -  HSSF  differences  later  in  this  report. 
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Figure  15.  Scatter  plot  of  instantaneous  observed  concentrations 
vs.  the  Mt.  Iron  multi-parameter  regression  equation  for  HSSF 
releases.  Data  have  been  adjusted  to  "95%  confidence"  and  the 
corresponding  LVDE  regression  equation  is  listed  at  the  top. 
Diagonal  lines  represent  values  within  a  factor  of  two  and  four  of 
agreement.  No  data  have  been  edited  when  calculating  these 
"goodness  of  fit"  grades. 
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Figure  16.    Same  as  figure  15,  except  data  are  time  averaged 
concentrations  from  HSSF  releases. 
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Figure  17.    Same  as  figure  15,  except  data  are  instantaneous 
concentrations  from  Lompoc  Valley  releases. 
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Figure  18.    Same  as  figure  15,  except  data  are  time  averaged 
concentrations  from  Lompoc  Valley  releases. 


5.   Conclusions  and  Recommended  Regression  Equations 

We  observe  significant  differences  between  the  LVDE  data  and 
the  Mt.  Iron  equations,  mainly  in  the  low  dose  range  where  the  Mt. 
Iron  values  are  lower  than  those  observed  in  LVDE.  We  speculate 
but  cannot  confirm  that  these  differences  may  stem  from  losses 
of  the  zinc  sulfide  aerosol  used  in  Mt.  Iron  due  to  gravitationally 
settling  and  impact  with  uneven  terrain.  Significant  improvement 
in  the  regression  predictions  are  not  obtained  by  adding  turbulence 
or  wind  speed  to  the  set  of  regression  parameters.  The  short  (1 
minute)  time  of  averaging  may  be  responsible  for  the  low 
correlation  with  turbulence.  Correlation  of  wind  speed  with  other 
regression  parameters  (e.g. AT)  most  likely  negates  improvement  in 
the  regression's  concentration  predictions.  Regression  predictions 
for  the  HSSF  release  site  are  significantly  different  from  those 
for  Lompoc  Valley  releases. 

Because  of  these  ambiguities,  we  recommend  validation, 
verification  and  certification  of  a  more  physically  based  diffusion 
model  for  cold  spill  modeling  at  Vandenberg.  If  a  verified 
regression  algorithm  must  be  used  in  the  interim,  we  recommend  the 
equation  given  in  fig.  16  for  predictions  of  hourly  averaged 
concentrations  resulting  from  continuous  releases  at  the  HSSF,  and 
the  equation  in  fig.  15  for  instantaneous  concentrations.  LVDE 
measurements  were  made  along  trajectories  between  the  HSSF  and 
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Lompoc  or  Miguelito  Canyon  and  only  during  daylight  hours. 
Therefore,  we  only  recommend  the  use  of  our  regression  results  for 
flows  in  this  general  direction  and  for  non-stable  boundary  layers. 
The  original  Mt .  Iron  equations  may  be  used  for  time  averaged 
concentrations  for  releases  from  SLC4  or  SLC6,  but  caution  should 
be  used  when  applying  them  to  long  ranges  (i.e.  off  base  locations) 
where  measurements  were  not  made.  Regressions  of  these  types 
should  not  be  used  for  atmospheric  conditions  or  locations  not 
included  in  the  measurement  domain.  Normally,  data  from  the  met 
tower  closest  to  the  release  site  should  be  used,  but  we  have  shown 
that  this  may  not  always  give  the  best  predictions,  due  to  passage 
of  the  plume  from  cloud  covered  to  clear  skies. 


6.   Basic  Equations 

Our  next  task  is  to  reduce  the  number  of  free  parameters  by 
using  the  non-dimensional  quantities  generally  accepted  in  the 
modern  literature.  We  must  first  introduce  a  few  fundamental 
equations. 

In  our  case,  a  simple  Gaussian  plume  model  predicts  a  centerline 
ground-level  concentration  normalized  by  release  rate  as 
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-^  = ,  (9) 

where  av  and  az  are  the  lateral  and  vertical  plume  parameters  for  a 
given  downwind  distance.  An  idealized  Gaussian  crosswind  shape 
also  predicts  that  the  surface  level  cross-wind  integrated 
concentration  (CWIC)  at  a  given  downwind  distance  can  be  calculated 
from 


/  C  dy  =  yf2n  ayCm     .  (10) 


If  the  plume  is  allowed  to  thoroughly  mix  through  the  depth  of  the 
well-mixed  boundary  layer,  h,  a  simple  mass  balance  within  the 
plume  reguires  that 


_  Uh  I  C  dy=,  (11) 

y  "      Q  ' 


where  Cy  is  the  normalized  CWIC.    Under  these  same  well  mixed 
conditions,  mass  balancing  the  Gaussian  plume  model  reguires  that 
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which,  in  turn,  forces  the  well-mixed  vertical  plume  parameter  to 
be 


o_ 


\ 


2  h     .  (13) 


While  there  is  ample  evidence  for  Gaussian  profiles  of  cross-wind 
concentration  from  continuous  surface  sources,  vertical  mass 
distributions  are  often  quite  non-Gaussian.  We  apply  this  equation 
only  as  a  qualitative  check  on  the  well-mixed  nature  of  the  plume. 


7.   The  Convective  Scaling  Approach 

The  convective  scaling  approach  has  been  shown  to  adequately 
predict  diffusion  from  surface  or  elevated  sources  within  buoyancy 
driven  boundary  layers  in  numerical  models  (Lamb,  1982) ,  tank  tests 
(Willis  and  Deardorff,  1981),  and  field  experiments  (Hanna,  1986). 
Most  applications  have  been  to  homogeneous  terrain,  but  Sakiyama 
and  Davis  (1987)  showed  success  in  terrain  of  varying  roughness. 

All  applications  scale  results  against  the  dimensionless  convective 
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scaling  distance, 


x  -  ^T    ■  <14) 


where  w*,  the  convective  scaling  velocity,  is  defined  as 


w.   = 


'Mil)1"  (is) 


and  where  H  is  surface  heat  flux,  g  is  gravitational  acceleration, 
and  T,  Cp  ,  p  are  average  temperature  within  the  boundary  layer,  and 
specific  heat  and  density  near  the  surface.  Note  that  X  can  be 
interpreted  as  the  ratio  of  the  plume  travel  time  (x/U)  to  an 
integral  time  scale  of  the  boundary  layer  turbulence  (h/w*) . 
Therefore,  convective  scaling  is  most  applicable  when  X  is  fairly 
close  to  unity.  As  X  strays  from  unity,  we  must  have  increasing 
faith  that  convective  turbulence  dominates  other  turbulence  sources 
at  scales  which  may  be  more  appropriate  to  the  plume  travel  time 
(e.g.  building  wakes  at  short  X,  mesoscale  wind  flow  variations  at 
large  X) . 

As  a  conservative  "lower  limit"  for  lateral  dispersion,  Briggs 
(1985)  suggests  the  following  non-dimensional  parameterization; 
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(1+2*) 


1/2 


(16) 


which  closely  approximates  Lamb's  numerical  model  and  Deardorff s 
tank  tests.  Similar  equations  are  presented  in  the  literature  for 
oz.  Other  equations  predict  the  normalized  CWIC,  Cy.  Both 
quantities  approach  limits  when  the  plume  becomes  well  mixed 
through  the  depth  of  the  boundary  layer  (roughly  X  >  1)  . 

The  basic  characteristics  of  eq.  16  adhere  to  the  theoretical 
limits  of  lateral  plume  growth  as  predicted  by  statistics;  Oy  ~  x 
for  short  distances  and  o  ~  x  '  at  long  range.  For  example, 
Draxler's  (1976)  classical  statistical  analysis  determines 
dispersion  parameters  from 
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(17) 


where  t  is  plume  travel  time,  T^  is  an  integral  time  scale 
proportional  to  Tl,  the  Lagrangian  time  scale.  Again,  the  plume 
grows  linearly  for  short  travel  time  and  as  t1/2  far  from  the 
source. 
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Field  experiments  such  as  Prairie  Grass  have  shown  that  eq.  16 
should  be  considered  as  a  lower  limit  to  lateral  diffusion  because 
low  frequency  turbulence  tends  to  increase  plume  spread.  Indeed, 
using  hundreds  of  hours  of  data,  Hanna  (1986)  showed  that  the 
following  simple  linear  equation  approximates  diffusion  of  elevated 
plumes  to  distances  of  50  km; 


-^  =  0.6*  .  (18) 

h 


A  surface  plume  release  would  tend  to  disperse  in  a  different 
fashion  for  short  X,  but  would  also  follow  a  linear  relationship  at 
distances  significantly  greater  than  unity.  We  must  emphasize  that 
this  equation  holds  on  average,  and  other  more  'controlled*  field 
experiments  show  growth  less  than  linear  for  large  X.  We  therefore 
consider  Hanna ' s  equation  to  be  an  upper  limit  for  lateral 
diffusion. 

For  X  <  4,  the  effects  of  'top  down  -  bottom  up'  diffusion  have 
been  observed  numerically  in  tank  tests  and  in  the  field. 
Materials  released  near  the  surface  tend  to  rise  or  'loft',  while 
the  height  of  maximum  concentration  of  materials  released  near  the 
top  of  the  boundary  layer  tends  to  lower.  Consequently,  Cywill  be 
less  than  unity  for  1  <  X  <  4,  and  surface  maximum  concentration 
will  be  less  than  that  predicted  with  eqs.  9  or  12.  The  reason  is 
given  by  the  structure  of  vertical  turbulence   in  convective 
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boundary  layers;  downdrafts  are  strong  and  spatially  confined, 
while  updrafts  are  weak  but  distributed  over  broad  regions.  As  the 
materials  become  well  mixed  in  the  vertical,  the  influence  of  this 
directional  asymmetry  wanes. 


7.   Variation  of  Convective  Scaling  Variables 

The  discussion  in  the  previous  section  assumes  horizontally 
homogeneous  conditions.  When  we  try  to  select  appropriate  w* ,  U, 
and  h  values  to  calculate  the  non-dimensional  convective  scaling 
distance,  X,  we  quickly  discover  significant  variations  in  those 
parameters  across  the  LVDE  measurement  domain.  Commonly,  LVDE 
releases  were  made  under  clouds.  However,  the  plume  pointed  inland 
and  measurements  were  usually  made  in  clear  skies.  In  Skupniewicz, 
et  al.  (1990,  1991b)  we  measured  and  modeled  changes  in  boundary 
layer  parameters  across  the  cloud  edge.  We  found  that  boundary 
layer  mean  windspeed  were  as  much  as  doubled  in  the  vicinity  of  the 
edge,  and  termed  this  accelerated  flow  a  "cloud  breeze".  We  also 
found  that  the  boundary  layer  height  increased  rapidly  within  a  few 
kilometers  of  the  edge,  but  the  rise  was  quickly  damped  and  the 
total  fractional  rise  was  small  due  to  the  thick  subsidence 
inversion.  In  most  cases,  that  rise  is  within  the  natural 
variability  of  the  boundary  layer  depth,  measured  at  different 
locations  or  times.  To  no  surprise,  measured  surface  heat  flux  was 
two  to  four  times  as  large  in  the  sun.    Using  this  heat  flux 
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change,  an  idealized  numerical  model  adequately  simulated  the  wind 
speed  and  boundary  layer  height  modifications  at  the  cloud  edge. 

Skupniewicz  et  al.  (1991a)  describes  all  meteorological  data 
available  from  LVDE.  Three  SODARs  were  located  along  the  plume 
trajectory;  a  permanent  facility  1  km  downwind  of  the  release  site 
(heretofore  identified  as  "900"),  a  temporary  unit  in  central 
Lompoc  ("WTP"),  and  a  mobile  unit  ("trailer").  During  most  of  the 
HSSF  releases  (the  first  5  of  8  measurement  days)  ,  the  mobile  SODAR 
was  positioned  along  the  trajectory  between  900  and  WTP.  During 
Lompoc  Valley  releases  (the  last  3  days) ,  the  trailer  was  at  the 
release  site,  and  there  was  no  intermediate  SODAR.  A  set  of  solar 
radiometers  was  located  at  the  WTP  site  and  near  the  900  facility 
from  which  surface  heat  flux  was  estimated  (see  Appendix  A)  . 
Figure  19  shows  the  fixed  SODAR  locations  with  average  plume 
trajectories  for  each  release  location. 

We  cannot  produce  a  detailed  picture  of  horizontal  variations 
in  w* ,  U  or  h  along  the  plume  trajectory  from  the  available 
measurements,  but  we  can  show  that  the  general  characteristics 
described  in  Skupniewicz  et  al.  (1991b)  were  present  during  LVDE. 
Figure  20  plots  w*  at  900  against  w*  at  WTP  for  the  times  of  LVDE 
tracer  measurements.  900  was  mainly  cloud  covered  while  WTP  was 
almost  always  sunny.  The  exceptions  are  16  Aug,  when  clouds  were 
present  at  Lompoc  during  some  of  the  tracer  measurements,  and  17 
Aug,  when  conditions  were  sunny  at  all  locations.     For  other 

47 


times,  we  see  a  twofold  increase  in  heat  flux  in  the  sun  which 
produces  a  30%  increase  in  w*  .  This  change  in  heat  flux  is  somewhat 
lower  than  that  observed  by  Skupniewicz  et  al.  (1991b),  and  we 
attribute  this  to  a  slightly  lower  boundary  layer  and  less  dense 
stratus. 

During  10-12  Aug,  all  three  SODARs  were  operating,  and  the 
mobile  one  was  located  in  clear  skies  a  few  kilometers  east  of  the 
cloud  edge.  Figure  21  shows  boundary  layer  height,  h,  measured 
concurrently  at  the  three  SODARs  during  LVDE  measurements  plotted 
against  their  average.  It  is  apparent  that  h  rises  significantly 
as  the  plume  exits  the  clouds,  but  lowers  to  roughly  its  original 
height  at  Lompoc .  The  lowering  is  contrary  to  the  numerical 
modeling  of  Skupniewicz  et  al.  (1991b)  based  solely  on  heating 
changes  between  cloudy  and  clear  skies,  but  consistent  with  some  of 
the  measurements  at  distances  well  inland  from  the  cloud  edge. 
SODAR  measurements  at  the  other  permanent  Vandenberg  facilities, 
located  at  the  coast  to  the  north  and  south  of  900  (mostly  cloud 
covered),  generally  agree  with  the  900  and  WTP  values.  While  high, 
the  trailer  values  are  within  the  range  of  variability. 
Unfortunately,  we  do  not  have  other  boundary  layer  heights  at  more 
inland  sites  to  better  define  h  in  the  clear  skies.  We  conclude 
that  while  the  h  increase  at  the  cloud  edge  is  significant,  the 
available  measurements  indicate  that  h  relaxes  to  its  original 
height  at  Lompoc,  the  eastern  edge  of  the  LVDE  domain. 
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Figure  19.  Topographic  map  of  the  LVDE  domain.  Lines  are  average 
plume  trajectories  for  each  release  site.  Permanent  SODARs  were 
located  at  Bldg.  900  and  WTP.  A  third  temporary  SODAR  was  located 
at  various  intermediate  positions  along  the  plume  trajectory. 
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Figure  20.  Surface  heat  flux  at  the  Water  Treatment  Plant  (+)  and 
the  release  site  (X) .  w.  at  the  Water  Treatment  Plant  (squares) 
and  the  release  site  (diamonds)  .  Data  are  one  hour  averages 
centered  on  the  time  of  release  for  each  LVDE  cross  wind  plume 
transect. 
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Figure  21.  Boundary  layer  heights  near  the  release  site  (Building 
900) ,  a  midpoint  along  the  plume  trajectory  (trailer) ,  and  in  the 
city  of  Lompoc  (WTP)  for  times  corresponding  to  plume  transects. 
Data  at  the  three  locations  are  plotted  against  their  average. 
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Figure  22.  Boundary  layer  vector  average  wind  speed  near  the 
release  site  (Building  900) ,  a  midpoint  along  the  plume  trajectory 
(trailer) ,  and  in  the  city  of  Lompoc  (WTP)  for  times  corresponding 
to  plume  transects.  Data  at  the  three  locations  are  plotted 
against  their  average. 
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Figure  23.  Boundary  layer  vector  average  wind  speed  near  the 
release  site  (Building  900),  a  midpoint  along  the  plume  trajectory 
(trailer),  and  in  the  city  of  Lompoc  (WTP)  for  times  corresponding 
to  plume  transects.  Data  at  the  three  locations  are  plotted 
against  the  vector  average  of  12  ft  wind  speeds  measured  at  24 
towers  within  20  km  of  the  release  site  multiplied  by  the  factor 
1.5. 


In  a  similar  fashion,  wind  speeds  averaged  through  the  depth 
of  the  boundary  layer  are  plotted  against  their  average  in  fig.  22. 
In  most  cases,  speeds  at  the  trailer  are  significantly  higher  than 
speeds  at  900  and  WTP,  which  generally  agree  with  each  other. 
Again,  this  pattern  agrees  with  the  conceptual  model  of  the  "cloud 
breeze"  proposed  by  Skupniewicz  et  al .  (1991b).  As  with  the 
boundary  layer  heights,  our  data  show  that  the  speed  enhancement 
was  negligible  at  the  eastern  edge  of  the  measurement  field. 

Wind  estimates  from  one  or  more  SODARs  were  often  unavailable 
due  to  the  degrading  effects  of  stratus  or  equipment  failure. 
Therefore,  we  would  like  to  use  tower  measurements  to  estimate 
boundary  layer  average  wind  speed.  Figure  23  plots  the  average  of 
all  operating  12  ft  wind  speeds  (maximum  24  towers)  against  the 
three  SODAR  measurements  discussed  above.  After  adjusting  the  12 
ft  speeds  by  the  factor  1.5,  the  characteristics  of  fig.  22  are 
duplicated  with  similar  scatter.  Such  an  ad  hoc  wind  speed 
estimate  is  advised  only  when  true  field  measurements  are  available 
for  calibration. 


8.   Convective  Scaling  Applied  to  LVDE 

Clearly,  the  parameters  that  determine  the  non-dimensional 
convective  scaling  distance,  X,  are  changing  along  the  trajectory 
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of  the  plume.  To  demonstrate  the  general  characteristics  of  plume 
parameters,  we  assume  horizontally  homogeneous  conditions  for  each 
plume  measurement.  Later,  we  allow  a  change  in  scaling  based  on 
the  position  of  the  cloud  edge  along  the  plume  trajectory. 

We  consider  only  hourly  averaged  plume  parameters  for  this 
analysis.  For  each  data  point,  time  averages  of  meteorological 
quantities  are  centered  upon  the  estimated  time  of  release,  after 
adjustment  by  the  approximate  plume  travel  time.  In  the  following, 
we  calculate  X  with  the  average  w*  measured  at  SODARs  900  and  WTP 
(fig.  20),  h  estimated  from  the  average  of  the  three  SODARs  (fig. 
21),  U  estimated  from  the  adjusted  base  average  wind  speed  (fig. 
23)  . 

Figure  24  shows  growth  in  the  lateral  plume  dimension,  Oy. 
Also  shown  are  the  Briggs  (1985)  and  Hanna  (1986)  formulations  we 
consider  to  be  the  lower  and  upper  limits.  We  see  that  most  data 
lie  between  the  two  limits.  The  few  data  points  below  X  =  1  agree 
with  the  lower  limit.  The  plume  grows  rapidly  near  X  =  1,  then 
slowly  to  a  distance  of  roughly  X  =  5.  The  plume  again  expands 
rapidly  beyond  X  =  6. 

For  each  transect  ay  was  calculated  directly  with  the  "moment 
method"  and  indirectly  by  applying  eq.  10  ("maximum  method"). 
Averaging  was  performed  with  eqs.  7-8.  The  pattern  described  above 
is  repeated  for  both  methods.   Since  very  little  difference  is 
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observed,  we  use  the  more  direct  moment  method  av  from  this  point 
forward.  Data  obtained  during  the  Lompoc  Valley  release  are 
depicted  as  shaded  symbols.  We  found  differences  between  the 
release  locations  negligible,  so  we  keep  all  data  grouped  together. 

Figure  25  shows  the  normalized  CWIC,  Cy,  calculated  directly 
with  eq.  11,  and  indirectly  by  applying  eq.  10  to  eq.  11.  Data 
below  X  =  6  tend  to  cluster  about  unity,  implying  that  the  plume  is 
is  well  mixed  in  the  vertical.  Plume  lofting  would  be  indicated  by 
a  decrease  in  Cy  f  rom  X  =  1  to  X  =  4 ,  and  we  see  no  evidence  of  this 
effect.  Unexpectedly,  Cydata  beyond  X  =  6  are  significantly  larger 
than  unity.  This  implies  an  increase  in  plume  mass  measured  at  the 
surface,  contrary  to  intuition.  We  will  offer  a  plausible 
explanation  for  this  anomaly  later. 

We  identify  Lompoc  Valley  releases  in  fig.  25  with  shaded 
symbols,  and  we  see  no  reason  to  separate  Cydata  into  subsets  based 
on  their  respective  release  locations. 

To  further  test  the  assumption  of  a  vertically  well  mixed 
plume,  we  estimate  az,  the  vertical  plume  dimension,  with  eq.  9  and 
take  the  ratio  of  oz  and  the  theoretical  limit  calculated  from  eq. 
13.  Figure  25  shows  values  near  or  above  unity  with  no  upward 
trend,  indicating  that  the  plume  has  arrived  at  its  limiting 
vertical  dimension,  if  no  plume  lofting  has  occurred.  Again  we  see 
a  change  at  roughly  X  =  6,  where  the  plume  appears  to  be  shorter 
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than  its  limiting  value.   This  trend  is  undoubtedly  an  artifact  of 
the  increased  CWIC  measured  at  the  surface  and  noted  above. 

Figure  26  shows  maximum  concentrations  normalized  by  release 
rate  calculated  with  three  different  methods.  In  contrast  to  the 
ay  methods,  the  "maximum  method"  is  the  more  direct  measurement 
while  the  "moment  method"  implies  maximum  concentration  from  eq. 
10.  A  third  method  shown  in  fig.  26  uses  eq.  12,  assuming  the 
plume  is  well  mixed  in  the  vertical.  The  maximum  and  moment  method 
closely  agree,  while  the  "well  mixed  method"  agrees  with  the  others 
only  for  distances  less  than  X  =  6.  Beyond  that  point  the  "well 
mixed  method"  is  lower  than  the  others,  presumably  due  to  the  rapid 
lateral  expansion  of  the  plume.  We  use  the  direct  measure  of 
maximum  concentration  from  this  point  forth. 

The  Lompoc  Valley  releases,  indicated  in  fig.  26  with  shaded 
symbols,  are  somewhat  lower  and  tend  to  segregate  from  the  HSSF 
releases.  We  recall  observing  lower  Lompoc  Valley  centerline 
concentrations  in  the  regression  analysis.  Due  to  the  large 
scatter,  data  sparsity,  and  lack  of  similar  differences  in  ay  or  Cy, 
we  consider  the  lower  Lompoc  Valley  Cq/Q  values  a  statistical 
anomaly  and  thus  will  keep  all  measurements  grouped  together. 
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Figure  24.  One  hour  averaged  lateral  plume  spread  as  a  function  of 
average  nondimensional  distance.  a  is  estimated  from  moment 
calculations  of  the  cross  wind  mass  distribution  and  a  is 
estimated  from  eg.  10.  For  each  case,  upper  line  is  the  upper 
limit  as  defined  by  eg.  18  and  the  lower  line  is  the  lower  limit, 


eg.  16.   For  o 


ym ' 


shaded  symbols  indicate  Lompoc  Valley  releases. 
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Figure  25.  Normalized  cross  wind  integrated  concentration,  C  , 
calculated  directly  (eq.  11) ,  and  estimated  by  substitution  of  eq. 
10  into  eq.  11  (C  )  as  a  function  of  average  nondimensional 
distance.  For  C  ,  shaded  symbols  indicate  Lompoc  Valley  releases. 
Also,  ratio  of  eq.  9  estimate  of  vertical  plume  spread,  o  ,  to  its 
well-mixed  limit  (eq.  13) . 
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Figure  26.  Normalized  surface  maximum  concentration  measured 
directly,  C0X/Q,  and  approximated  from  eq.  10  (C^Q)  as  a  function 
of  average  nondimensional  distance.  For  C^Q,  shaded  symbols 
indicate  Lompoc  Valley  releases.  Also,  Coe/Q  is  estimated  assuming 
plume  is  well-mixed  through  the  boundary  layer  (eq.  12)  . 


9.   Two  Zone  Convective  Scaling 

The  approach  we  take  is  to  consider  the  domain  as  two  zones 
demarcated  by  the  cloud  edge.  For  each  measurement,  we  have 
calculated  the  approximate  location  of  the  cloud  edge  along  the 
plume  trajectory  based  on  hourly  GOES  satellite  images.  Appendix 
B  details  our  methodology  for  this  procedure.  We  have  already 
described  changes  in  w* ,  U,  and  h  attributed  to  the  cloud  edge. 
Here,  we  only  consider  a  step  change  in  w* at  the  edge,  maintaining 
the  average  values  of  U  and  h  used  previously.  We  would  prefer  to 
allow  all  scaling  parameters  to  vary,  but  we  cannot  detail  the 
changes  in  U  or  h  along  the  plume  trajectory  based  on  our 
measurements.  Such  a  model  input  requirement  would  certainly  be 
difficult  from  an  operational  standpoint. 

Refer  to  fig.  27.  Let  curves  1  and  2  define  plume  growth 
referenced  to  non-dimensional  distance  calculated  from  scaling 
parameters  under  the  stratus  (zone  1)  and  in  the  sun  (zone  2), 
respectively.  Curve  3  defines  a  plume  which  crosses  the  zone  1-2 
boundary.  In  our  case,  this  boundary  is  the  cloud  edge,  Xc.  We 
assume  that  the  plume  grows  with  zone  2  scaling,  as  if  it  had  been 
released  from  a  virtual  source,  Xv,  some  distance  between  X=0  and 

Xq  . 
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Figure  27.  Schematic  of  two  zone  convective  scaling  for  lateral 
plume  spread.  For  LVDE,  zone  1  is  cloud  covered  and  zone  2  is 
cloud  free.  Curve  1  defines  plume  growth  in  zone  1  coordinates, 
curve  2  defines  growth  in  zone  2  coordinates.  Two  zone 
parameterization  follows  curve  1  in  zone  1  and  curve  3  in  zone  2. 


X  is  the  cloud  edge  and 


is  the  virtual  source  (e.g.  eq.  23) . 


The  function  defining  curve  3  can  take  many  forms,  depending 
on  the  choices  for  curves  1  and  2  .  One  choice  is  the  Briggs 
formulation  (eq.  16) .  Calculating  X  with  scaling  parameters  valid 
in  zone  1,  curve  2  is  defined  by 


h  (1+2FX) 


1/2 


(19) 


where  F  is  the  ratio  of  the  respective  scaling  parameters, 


w. 


F  =    — 2—^-      ,  (20) 


U.h, 


In  our  case,  we  assume  U  and  h  are  constant  (averaged  over  the 
domain)  ,  so  F  =  w.2/w.,.   Curve  3  is  defined  by  the  set  of  equations 


X  <.  X  (21) 


h  {1+2X)1/2 


ov7      0.6F(X-XV) 
-/I    =   1 x>Xr  (22) 

h  (1+2F(X-XJ1/2 
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Xv   =      iF   1]  Xc     ,  (23) 


The  equation  for  Xv  is  obtained  by  matching  the  two  previous 
equations  at  Xc.  One  can  see  that  F  can  take  any  value,  so  we  can 
equivalently  choose  to  model  a  plume  entering  a  stratus  covered 
domain  from  the  clear  skies. 

We  have  stated  earlier  that  eq.  16  should  be  considered  a 
lower  limit.  We  suggest  a  generic  parameterization  which  fits  most 
lateral  plume  growth  formulations  suggested  in  the  literature: 


-^  =  AXa  X<Xn,  (24) 


^=AX01la-b)Xb  XQ1<   X  ±  Xc  (25) 


^  =  AXQ1{a-b)  (F(X-XV)  )  b       X>XC  (26) 
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Xv  =  X-  — 

v        F 


(  x^\ 

\   Xq2  ) 


(^) 


(27) 


where  A, a, b, Xqi ,and  XQ2are  free  parameters  and  Xc >  Xqi  Most  short 
range  experiments  and  numerical  modeling  results  would  reguire  A  = 
0.6  and  a  =  1  (i.e.,  Prairie  Grass).  Since  we  have  few  data  for  X 
<  1,  we  use  these  values.  b  determines  growth  at  large  distances, 
and  is  in  the  range  of  0.5  to  1.  We  prefer  to  use  b  =  2/3  as 
suggest  by  Deardorff  and  Willis  (1975),  Lamb  (1979),  and  Briggs 
(1985)  . 

Xoiand  XQ2are  the  distances  at  which  plume  growth  transitions 
from  Xa  to  X  in  zones  1  and  2,  respectively.  Conceptually,  the 
transition  occurs  when  the  width  of  the  plume  has  exceeded  the 
dominant  scale  of  the  convective  energy  producing  eddies. 
Selecting  a  transition  point  value  of  0.18  adequately  reproduces 
the  numerical  modeling  results  of  Lamb  (1979)  and  the  tank  tests  of 
Deardorff  and  Willis  (1975).  Briggs  (1985)  suggests  that  a 
transition  at  X  =  0.6  provides  a  better  match  to  field  experiments. 

We  contend  that  the  transition  points  of  zone  1  and  zone  2  can 
be  different.  Pasquill  and  Smith  (1983)  shows  that  we  can 
calculate  plume  spread  by  integrating  the  Lagrangian  turbulence 
spectrum  after  applying  a  low  pass  filter  at  a  frequency  inversely 
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proportional  to  the  plume  travel  time.  Therefore,  the  transition 
point  should  occur  when  the  travel  time  is  near  the  Lagrangian  time 
scale.  Thinking  in  terms  of  spectral  similarity  in  the  convective 
boundary  layer  (e.g.  Hojstrup,  1982),  the  peak  in  the  Eulerian 
lateral  energy  spectrum  occurs  near 

nm  =   -*—    =  — ~    =0.5   ,  (28) 


th 


U  TmU 


where  fm  is  the  frequency  corresponding  to  a  dominant  eddy  time 
scale  Tm.   By  definition, 


X    -     KlX    -     Kltl  (29) 


where  t]_  is  plume  travel  time  to  X]_.  So,  it  follows  that 


Xl   -  Z-Zk  «  il    .  (30) 


Since  we  use  an  average  wind  speed  and  w*iis  smaller  than  w*2  we 
contend  X^  should  be  smaller  than  X2  •  The  proportionality  constant 
in  eq.  30  depends  on  the  Eulerian-Lagrangian  time  scale  ratio  which 
is  difficult  to  measure  or  estimate,  so  without  detailed  spectral 
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information  we  state  only  that  the  constant  may  also  be  different 
in  zones  1  and  2.  Fortunately,  the  model  is  not  extremely 
sensitive  to  these  parameters.  Using  reasonable  values  obtained 
from  the  literature,  we  choose  X^ =  0.18  (approximating  Lamb,  1979) 
and  X2  =  0.60  (Briggs,  1985). 

With  these  assumptions,  the  generalized  equations  reduce  to 


-^  =  0.6X         X   <  0.  18  (31) 


-^  =0.34  X2n  0.18  <  X  ±  Xr  (32) 

h 


-^  =  0.34  (F(X-XV)  )2/3  X  >   Xc  (33) 


Xv  =      (F-°-55)  Xc       ,  (34) 

F 


10.   Convective  Scaling  Model  Comparison  -  General 


In  this  section,  we  compare  measurements  of  ay,  Cq/Q,  and  Cy 
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against  the  two  zone  model  proposed  above  and  four  models  outlined 
earlier.  For  all  cases,  Cq/Q  and  Cy  are  calculated  from  eqs.  10-12 
after  calculation  of  ay.  The  four  test  models  are  1)  Briggs'  eq. 
16  obtained  from  scaling  parameters  in  zone  1  and  thus  considered 
the  "lower  limit",  2)  Hanna ' s  eq.  18  obtained  from  scaling 
parameters  in  zone  2,  considered  the  "upper  limit",  3)  Briggs1 
(1985)  "best  fit"  recommendation  (eqs.  24  and  25  with  A  =  0.6,  a  = 

1,  b  =  2/3,  and  X]_  =  0.6)  obtained  from  scaling  parameters  in  zone 

2,  considered  an  intermediate  choice,  and  4)  Draxler's  (1976) 
statistical  approach  (eq.  17)  applied  to  LVDE  data. 

We  have  only  briefly  described  the  statistical  approach. 
This  method  requires  good  measurements  of  turbulence  intensity,  oq. 
We  use  the  10  minute  average  value  measured  by  Vandenberg's  3-axis 
Gill  anemometer  at  Building  900,  the  bivane  mounted  on  Tower  057, 
or  the  trailer's  sonic  anemometer,  depending  on  the  release  site 
and  data  availability.  The  method  is  only  valid  for  plume  travel 
times  within  roughly  an  order  of  magnitude  of  the  Lagrangian  time 
scale.  This  method  gives  reasonable  results  only  after  excluding 
the  above  described  long  range  data,  where  rapid  cloud  expansion 
occurs.  Figure  28  shows  the  results  after  this  editing.  The 
derived  integral  time  scale,  ti=  330  sec,  is  within  the  range  of 
values  calculated  by  Draxler  for  surface  releases. 
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Figure  28.  Draxler's  (1976)  statistical  approach  applied  to  LVDE. 
+  are  measured  values  and  squares,  line  are  eq.  17  with  t.  =  33.0 
seconds.  Due  to  rapid  cloud  growth  at  long  distances  from  the 
source,  only  data  with  travel  times  less  than  2000  seconds  were 
considered. 


Figures  29  and  30  show  the  comparison  for  <jy .  Models  1  (lower 
limit)  and  2  (upper  limit)  underpredict  and  overpredict  plume 
spread  in  a  similar  fashion  to  fig.  24,  where  all  scaling 
parameters  were  averaged.  Model  3  (intermediate)  slightly 
overpredicts  av  at  short  distances  (X  <  6)  ,  and  underpredicts  for 
X  >  6.  Models  4  (Draxler)  and  5  (two  zone)  do  a  good  job  at  short 
range,  then  underpredict  at  long  range.  The  deficit  at  long  range 
is  eliminated  for  the  two  zone  model  only  if  the  ay  predictions  are 
doubled.   We  show  this  arbitrary  adjustment  in  the  figures. 

We  might  conclude  that  either  model  4  or  5  may  be  used  for 
distances  less  than  X  =  6.  We  caution  use  of  model  4  for 
meteorological  conditions  or  locations  different  from  those  of 
LVDE.  This  semi-empirical  model  depends  critically  upon  the 
integral  time  scale  chosen.  There  is  no  guarantee  that  the  time 
scale  is  valid  for  other  conditions  or  locations.  oq  was  measured 
directly  in  the  path  of  the  plume  with  fast  response  instruments. 
Using  the  slower  response  tower  wind  vanes  at  other  locations  will 
reduce  oq  and  the  predicted  plume  dimensions.  On  the  other  hand, 
the  two  zone  model  is  purely  theoretical  and  does  not  depend  on 
direct  turbulence  measurements. 
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Figure  29.  Ratio  of  predicted  to  measured  lateral  plume  spread  as 
a  function  of  nondimensional  distance  calculated  from  scaling 
parameter  measured  at  the  source  (zone  1)  .  Model  1  is  assumed 
lower  limit  (eq.  16) ,  model  2  is  assumed  upper  limit  (eq.  18) ,  and 
model  5  is  two  zone  (eqs.  31-34).  For  model  5,  small  squares  at  X 
>  6  are  model  5  predictions  multiplied  by  a  factor  two. 
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Figure  30. 
recomme 


30.   Same  as  figure  29,  except  model  3  is  Briggs'  (1985) 
ndation  and  model  4  is  eg.  17  statistical  estimate. 


Figures  31  and  32  show  C0/Q  predictions .  The  "limit"  models 
(1  and  2)  poorly  predict  maximum  concentration.  Model  4  (Draxler) 
slightly  overpredicts  centerline  concentration.  Models  3  and  5  do 
a  reasonable  job  at  all  ranges,  except  for  a  few  points  at  the 
furthest  range.  The  rapid  expansion  of  the  plume  at  X  >  6, 
mentioned  earlier,  is  not  reflected  in  the  Cq/Q  measurements.  If 
this  were  true,  we  would  see  a  sharp  decrease  in  centerline 
concentration  beyond  X  =  6.  We  offer  a  plausible  explanation 
later. 

Figures  33  and  34  show  Cy  predictions .  An  agreement  between 
observation  and  model  indicates  conservation  of  mass  at  the 
surface.  Again,  the  lower  and  upper  limit  models  are  severely 
biased  at  all  ranges.  Only  Draxler's  model  (4),  the  two  zone  model 
(5)  adeguately  predicts  Cy  for  X  <  6.  No  model  gives  adequate 
prediction  for  all  ranges  because  observations  of  Cy  sharply 
increase  beyond  X  =  6.  As  the  figures  show,  arbitrarily  doubling 
Oy  for  the  two  zone  model  at  X  >  6  also  reasonably  adjusts  the  Cy 
prediction. 

In  summary,  the  two  zone  formulation  adequately  predicts  Oy, 
Cq/Q,  and  Cy  as  judged  by  a  comparison  with  the  four  models  based  on 
homogeneous  conditions.  No  model  properly  predicts  all  quantities 
for  X  >  6.  Only  by  arbitrary  adjustment  of  predictions  at  these 
ranges  do  we  simulate  Oy   and  Cy  data.   Oddly,  this  adjustment  is 
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unneeded  for  Cq/Q. 

Figure  35  illustrates  the  following  argument.  When  we  closely 
examined  the  raw  data  for  points  beyond  X  =  6,  we  saw  very  "flat" 
or  bimodal  instantaneous  cross-wind  concentration  distributions. 
Referring  to  fig.  19,  all  these  data  were  measured  at  downwind 
distances  beyond  LaSalle  Canyon.  We  contend  that  the  shear  created 
by  the  ventilating  effects  of  Sloan  and  Miguelito  Canyons  rapidly 
expands  the  lateral  dimensions  of  the  plume  at  the  surface.  One 
need  only  briefly  visit  Miguelito  Canyon  to  realize  that  the 
daytime  flow  at  the  surface  is  distinctly  up-canyon.  Closer  to  the 
coast,  near  surface  up-valley  flow  is  negligible,  such  as  in  the 
cooler  Lompoc  Canyon.  The  low  level  northerly  surface  winds  of  the 
inland  side  canyons  are  in  sharp  contrast  to  the  general  west- 
northwest  flow  over  the  region,  as  indicated  by  plume  trajectories 
and  SODAR  winds.  Because  the  up-canyon  flow  is  so  shallow,  mass 
removed  from  the  "main"  plume  is  quickly  replaced  by  mass  from 
aloft  through  vertical  mixing.  The  net  result  is  a  wide,  flat 
surface  concentration  distribution.  The  cross-wind  integrated 
concentration  increases,  as  if  the  plume  were  gaining  mass. 

The  above  arguments  are  speculative,  but  are  consistent  with 
our  observations.  Only  a  mesoscale  numerical  model  which 
accurately  captures  the  main  features  of  these  up-canyon  flows 
could  unambiguously  determine  the  actual  causes. 
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Figure  31.  Ratio  of  predicted  to  measured  surface  maximum 
concentration  as  a  function  of  nondimensional  distance  calculated 
from  scaling  parameter  measured  at  the  source  (zone  1) .  Models  are 
identified  in  fig.  29. 
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Figure    32.       Same   as    fig.    31,    except    for   models    3,    4,    and    5. 
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Figure  33.  Ratio  of  predicted  to  measured  normalized  cross  wind 
integrated  concentration  as  a  function  of  nondimensional  distance 
calculated  from  scaling  parameter  measured  at  the  source  (zone  1) . 
Models  are  identified  in  fig.  29.  For  model  5,  small  sguares  at  X 
>  6  are  model  5  predictions  assuming  a  is  multiplied  by  a  factor 
two. 
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Figure  34.   Same  as  fig.  33,  except  for  models  3,  4,  and  5. 


Figure  35.  Schematic  of  topographically  enhanced  plume  spread. 
Near  surface  wind  shear  is  caused  by  ventilating  effects  of  side 
canyons  which  intersect  Lompoc  Valley.  Rapid  vertical  mixing 
maintains  maximum  concentration  while  enhancing  surface  level  cross 
wind  integrated  concentration. 


11.   Convective  Scaling  Model  Comparison  -  Statistics 

Statistics  were  calculated  for  the  models  discussed  above. 
Table  1  summarizes  the  results.   All  data  were  include  in  the  Cq/Q 
statistics.   Data  at  X  >  6  were  omitted  from  ay  and  Cv  statistics. 
The  general  remarks  of  the  previous  section  are  quantified  in  these 
results.   A  few  additional  comments  follow. 
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Table  1. 

Stati 

Statistic 

Mod 

°Y 

Lowe 

o/pi  average 

0.72 

o/p  median 

0.69 

FB 

0.44 

GFB 

0.59 

q.  range2 

0.44 

factor  2 

0.74 

factor  4 

1.00 

NMSE 

0.46 

c0/Q 


o/p  average 

2.40 

o/p  median 

1.80 

FB 

-.37 

GFB 

-.02 

q.  range 

1.  39 

factor  2 

0.62 

factor  4 

0.87 

NMSE 

0.51 

CY 

o/p  average 

1.73 

o/p  median 

0.70 

FB 

0.42 

GFB 

0.28 

q.  range 

0.60 

factor  2 

0.74 

factor  4 

0.90 

NMSE 

0.74 

Draxler 

Two  Zone 

1.17 

0.99 

0.93 

0.98 

-.03 

0.10 

-.19 

0.29 

0.74 

0.67 

0.90 

0.84 

1.00 

1.00 

0.21 

0.29 

Statistical  Convective  Scaling  Model  Comparison 

Mod  2  Mod  3      Mod  4      Mod  5 

Upper  Med. 

1.80  1.33 

1.58  1.33 

-.50  -.17 

-2.4  -3.5 

0.98  0.68 

0.65  0.87 

0.97  1.00 

0.64  0.21 


0.89  1.28  1.37  1.77 

0.55  1.02  1.22  1.19 

0.55  0.19  -.15  -.16 

0.02  0.00  -.01  -.01 

0.50  1.05  1.07  0.96 

0.54  0.68  0.66  0.76 

0.85  0.97  0.97  0.92 

1. 15  0.48  0.43  0.43 


4.27  3.14  2.48  2.44 

1.57  1.29  1.05  0.95 

-.45  -.21  0.00  0.12 

****  -1.3  -.43  -.07 

2.25  1.13  0.98  1.12 

0.61  0.65  0.74  0.71 

0.90  0.94  0.97  0.97 

0. 69  0.44  0.35  0.49 


1  o/p  is  the  ratio  of  observed  data  to  prediction 

2  q.  range  is  the  interquartile  range  between  25  and  75%  of  median 
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Since  these  data  are  geometric  (ratios) ,  the  median  value  is 
more  appropriate  than  the  average.  Fractional  bias  (FB)  is 
negatively  correlated  with  the  median.  Model  overpredictions 
(median  >  1)  are  usually  associated  with  negative  FB.  We  see  that 
models  4  and  5  are  clearly  superior  in  terms  of  median  and  FB, 
meaning  they  most  accurately  reproduce  the  ensemble.  Geometric 
fractional  bias  (GFB)  indicates  trends  in  the  bias.  One  can 
calculate  significant  FB  with  minimal  GFB  (e.g.  Cq  predictions)  or 
vice  versa.  We  see  large  negative  GFB  values  with  modest  FB  for 
model  3  ay  and  Cy  predictions.  This  indicates  that  many  more  data 
points  over-predict  than  underpredict ,  even  though  the  magnitude  of 
that  over-prediction  is  small. 

The  factor  grading  was  explained  earlier.  Interquartile  range 
(q.  range)  is  the  actual  range  of  values  corresponding  to  the 
predictions  falling  within  a  factor  of  two  of  the  observations.  Q. 
range  tends  to  correlate  with  the  bias,  and  should  only  be  used  as 
a  measure  of  the  relative  merit  of  two  models  with  equal  bias. 
Reviewing  the  table,  the  factor  grading  is  quite  insensitive  to  the 
model  quality.  As  discussed  earlier,  it  is  sensitive  to  the  number 
of  data  in  the  ensemble.  Therefore,  we  recommend  avoiding  factor 
grading  for  these  type  of  analyses. 

Normalized  mean  square  error  (NMSE)  is  a  good  measure  of  the 
scatter  of  data,  regardless  of  bias.   The  table  indicates  that 
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models  3,4  and  5  have  comparable  data  scatter  for  all  quantities. 
Ideally,  we  would  liked  to  have  seen  significantly  reduced  NMSE  for 
the  two  zone  model.  Unfortunately,  requiring  cloud  edge  location 
adds  degrees  of  freedom  to  the  model,  and  likewise,  scatter.  This 
most  likely  offsets  any  scatter  reduction  realized  with  the  added 
physics. 
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12.   Comparison  of  Regressive  and  Convective  Scaling  Results 

The  regression  analyses  of  sections  2-5  only  predict 
centerline  concentration.  Those  regressions  were  conducted 
seperately  for  the  HSSF  and  Lompoc  Valley  releases.  Table  2 
applies  the  Mt.  Iron  equation  and  the  recommended  time  averaged 
regression  equation  (see  section  5)  to  the  full  set  of  data.  The 
Two  Zone  model  results  are  repeated  from  table  1. 

We  see  comparable  statistics  for  the  LVDE  regression  and  Two 
Zone  models.  The  Mt.  Iron  regression  performs  very  poorly  in  terms 
of  scatter  (NMSE)  and  predicts  a  very  low  median  value  compared  to 
the  observations.  A  tendency  towards  unreasonably  low  predictions 
at  greater  range  was  discussed  in  the  regression  analysis.  Oddly, 
the  Mt.  Iron  fractional  is  negative,  contrary  to  the  lower  median 
predictions.   This  is  due  to  a  few  very  large  predictions. 

The  mutual  agreement  between  the  LVDE  regression  equation  and 
the  theoretically  based  two  zone  model  further  support  the 
contention  that  the  Mt.  Iron  equation  should  not  be  used  for 
releases  from  HSSF.  Taken  at  face  value,  we  may  be  tempted  to 
conclude  that  the  LVDE  regression  and  Two  Zone  model  are  equivalent 
predictors.  We  must  realize  that  the  LVDE  regression  is  being 
compared  to  the  very  same  data  from  which  it  was  derived.  The 
other  two  predictions  are  independent  of  the  data.  Given  another 
set  of  data  with  different  background  conditions,   the  LVDE 
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regression  could  perform  as  poorly  as  the  Mt .  Iron  regression  does 
here.  The  two  zone  model  is  only  confined  by  the  physics  from 
which  it  was  derived,  and  can  therefore  be  generalized  to  a  wide 
variety  of  conditions. 
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Table  2.   Convective  Scaling  vs.  Regression  Model  Comparison 


regression 

convective  scaling 

c0/Q 

Mt.  1 

ron 

LVDE 

Two  Zone 

o/p  average 

1.33 

1.50 

1.77 

o/p  median 

0.63 

1.07 

1.  19 

FB 

-.25 

-.04 

-.16 

GFB 

0.02 

-.01 

-.01 

q  range 

0.99 

0.87 

0.96 

factor  2 

0.  36 

0.75 

0.76 

factor  4 

0.68 

0.95 

0.92 

NMSE 

2  .43 

0.46 

0.43 
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13.  Conclusions  and  Recommended  Convective  Scaling  Equations  (Cq/Q) 

Application  of  convective  scaling  to  LVDE  reasonably  agrees 
with  the  observations  and  theories  of  other  investigators. 
Considering  the  partially  cloud  covered  LVDE  domain  as  two  zones 
with  different  w*  improves  the  predictions.  The  accuracy  of  this 
formulation  is  comparable  to  statistical  predictions  using  actual 
turbulence  measurements  and  integral  time  scales  deduced  from  the 
tracer  measurements.  While  convective  scaling  predictions  of 
centerline  maximum  concentration  are  satisfactory  for  all  ranges, 
predictions  of  lateral  plume  dimension  and  surface  level  cross-wind 
integrated  concentration  fail  for  X  >  6.  We  attribute  this  failure 
to  topographic  flow,  namely,  the  ventilating  effects  of  canyons 
crossing  the  mean  plume  trajectory.  We  speculate  that  centerline 
concentrations  are  not  reduced  because  the  low  level  nature  of  this 
up-canyon  flow  allows  for  rapid  mixing  of  mass  from  the  unaffected 
plume  aloft. 

We  recommend  using  the  two  zone  model  of  eqs .  31-34,  when  the 
downwind  position  of  the  cloud  edge  relative  the  source  is  known. 
This  information  can  be  obtained  from  satellite  imagery,  but  errors 
are  large  and  the  procedure  is  not  trivial.  We  recommend  visual 
estimates.  Observers  with  local  knowledge  of  the  terrain  and 
landmarks  should  be  able  to  locate  the  cloud  edge  with  more 
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accuracy  than  satellite  images.  If  funds  become  available,  we  also 
recommend  that  all  meteorological  towers  be  fitted  with  low  cost 
photokvoltaic  radiometers.  Such  a  network  would  give  an  adequate 
definition  of  the  cloud  edge  position.  On  the  other  hand,  we 
recommend  installation  of  an  accurate  set  of  short  and  long  wave 
radiometers  at  an  inland  met  tower  for  clear  sky  heat  flux 
estimates.  WT019,  WT014,  or  a  new  tower  situated  in  Lompoc  would 
be  good  candidates.  The  present  radiometers  at  WT301  should  be 
maintained  and  factory  calibrated  periodically.  These  will  supply 
adequate  heat  flux  estimates  under  the  clouds. 

In  lieu  of  cloud  edge  information,  we  recommend  using  the 
statistical  model  (eq.  17)  for  releases  near  the  Lompoc  Valley, 
during  "normal"  sea  breeze  conditions.  For  this  purposel ,  we 
suggest  using  the  10  minute  turbulence  measurements  from  Building 
900  (level  z  -  0  in  the  Building  900  DASS  data  file)  .  We  have 
already  discussed  the  problems  with  using  this  formula  at  other 
locations  or  in  different  meteorological  conditions. 

For  totally  or  partially  cloud  covered  cases  with  no 
information  on  the  cloud  edge  position,  we  recommend  using  model 
1  (eq.  16),  as  suggested  by  Briggs  (1985).  w*  should  be  calculated 
under  the  stratus.  This  will  give  a  conservative  prediction  of 
centerline  concentration.  Based  on  the  results  of  this  analysis, 
one  may  be  tempted  to  use  model  3  which  would  give  lower  centerline 
concentrations.   However,  we  believe  model  3  is  most  applicable  to 
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clear  sky  conditions.  The  reasonable  agreement  between  model  3 
and  the  data  was  mainly  due  to  the  fact  that  the  cloud  edge  was 
usually  close  to  the  source. 

For  totally  clear  skies,  model  5  (two  zone)  reduces  to  model 
3.   Therefore,  either  model  will  produce  the  same  results. 

All  of  these  models  require  boundary  layer  depth  and  wind 
speed  input.  We  showed  that  these  data  are  highly  variable  for 
LVDE  along  the  plume  trajectory.  Since  we  could  not  unambiguously 
define  a  function  of  this  variability,  we  used  averaged  values 
which  produce  reasonable  results.  This  approach  can  be  carried  too 
far,  however,  if  data  are  included  in  the  average  which  are  far 
removed  from  the  plume  trajectory.  For  boundary  layer  height,  we 
did  not  include  SLC6  or  Building  1764  data  in  the  average.  For 
wind  speed,  we  used  the  base  average  of  12  ft  winds,  but  calibrated 
that  average  against  actual  boundary  layer  averaged  wind  speeds 
from  the  SODARs  along  the  plume  trajectory.  We  do  not  recommend 
using  our  calibration  factor  for  operational  purposes.  Our 
recommendation  is  to  use  boundary  layer  height  and  average  wind 
speed  from  the  SODAR  closest  the  release  point.  For  HSSF  or  SLC4, 
that  SODAR  is  at  building  900.  For  north  base,  that  SODAR  is  at 
building  1764.  Boundary  layer  height  should  be  interpreted  from 
the  shadowgraphs  as  the  center  of  the  elevated  backscattered  power 
maximum  (BL2  in  Skupniewicz  et  al.  1990,  see  this  reference  for 
further  interpretative  procedures) .    For  wind  speed,  perform  a 
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simple  average  of  the  speeds  below  the  boundary  layer  height.  We 
used  a  vector  average,  but  this  procedure  is  cumbersome  and  can 
give  unreasonably  low  values  in  highly  sheared  conditions. 

These  procedures  should  be  well  suited  to  either  a  planning  or 
emergency  situation.  All  necessary  data  input  should  be  easily 
acquired  within  a  matter  of  minutes.  The  eguations  can  be 
programmed  into  any  portable  computer,  and  results  obtained  with 
negligible  run  time. 


14.   Recommended  Convective  Scaling  Equations  (ay,  Cy) 

If  the  lateral  plume  dimensions  or  total  ground  level  mass  are 
required,  one  must  consider  the  topography.  The  above 
recommendations  hold  for  releases  near  the  mouth  of  the  Lompoc 
Valley  during  "normal"  sea  breeze  conditions,  provided  that  the 
downwind  receptor  is  west  of  Sloans  Canyon.  At  further  distances, 
we  cannot  recommend  an  analytical  model.  When  an  operational 
numerical  model  becomes  available,  it  should  be  tested  against 
these  data  to  see  if  the  rapid  plume  expansion  and  enhanced  ground 
level  CWIC  are  adequately  simulated. 


90 


15.   Recommended  Instrumentation  Upgrades  and  Future  Work 

We  have  already  made  the  recommendation  of  adding  solar 
radiometers  to  the  meteorological  tower  network  to  define  the  cloud 
conditions  and  provide  radiation  data  for  calculating  surface  heat 
flux.  We  additionally  recommend  research  into  application  of 
doppler  RADAR  technology  at  Vandenberg  for  the  purpose  of  defining 
boundary  layer  heights  and  winds.  One  possible  candidates  is 
"classical"  doppler  weather  RADAR  (e.g.  NEXRAD) .  It  may  be 
possible  to  tap  the  data  from  one  of  the  Vandenberg  operational 
RADARs  to  produce  doppler  wind  estimates.  However,  wind  estimates 
from  these  low  frequency,  low  elevation  RADARs  are  subject  to  error 
from  ground  clutter,  and  the  technique  may  not  be  feasible  at 
Vandenberg.  High  frequency  vertical  profiling  RADARs  are  rapidly 
becoming  available.  We  recommend  a  model  that  not  only  profiles  in 
the  vertical,  but  also  scans  at  off-zenith  angles.  We  have 
emphasized  the  horizontal  variability  of  boundary  layer  heights  and 
winds  in  our  research,  and  we  feel  that  any  expensive  remote 
sensing  research  efforts  or  purchases  at  Vandenberg  should  address 
this  problem. 

We  recommend  continued  development  of  a  general  numerical 
diffusion  model  for  operational  application  at  Vandenberg  and  other 
Air  Force  facilities.  These  models  have  a  wide  range  of 
sophistication  which  highly  correlates  with  "run  time".    For 
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Vandenberg,  the  degree  of  sophistication  reguired  is  large. 
Therefore  we  recommend  parallel  development  of  simpler  models  to  be 
used  in  emergency  situations.  For  example,  parts  of  this  analysis 
can  be  applied  in  simple,  fast,  comprehensive  Air  Force  models  such 
as  AFTOX  or  ADAM.  Continued  development  of  such  models  is 
critical . 


Appendix  A.  Convective  Scaling  Velocity 

To  estimate  the  convective  scaling  velocity  under  clear  or 
cloudy  conditions  we  must  know  the  on-site  inversion  base  height, 
as  well  as  the  surface  heat  flux.  We  assume  that  the  inversion 
base  height  is  eguivalent  to  the  average  height  of  thermal 
penetration,  interpreted  as  the  maximum  elevated  backscatter  power 
return  from  the  SODARs.  Two  elevated  echo  layers  are  sometimes 
observed  for  stratus  cases,  below  and  above  the  clouds.  We  use 
the  center  of  the  thin  cloud  top  echo  layer  as  the  inversion  base 
height.  These  are  estimated  by  visually  interpreting  the  SODAR 
facsimile  records. 

We  also  developed  an  algorithm  to  estimate  surface  heat  flux 
from  the  solar  and  infrared  sensors  we  had  placed  at  these  two 
sites.  At  present  VBG  has  only  one  set  of  radiation  sensors. 
Thus,  we  also  included  a  routine  to  estimate  downwelling  solar  and 
thermal  radiation,  if  measured  data  is  not  available.  This  routine 
reguires  a  fractional  cloud  coverage  estimate  as  well  as  other 
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commonly  available  input  data,  such  as  inversion  height,  screen 
level  wind  speed,  temperature,  and  relative  humidity.  The 
remainder  of  this  section  describes  this  routine,  as  well  as  the 
heat  flux  algorithm. 

All  solar  radiation  models  require  current  date  and  time  to 
compute  sun  angle,  \p .  This  is  done  in  module,  "SOLALT",  which  also 
computes  sunrise/set  times.  We  further  modified  an  improvement  of 
the  simple  ASHRAE  model  for  downwelling  solar  direct  and  diffuse 
radiation  (Iqbal,  1983)  by  supplying  sine  curve  fits  to  seasonal 
adjustments  for  apparent  extra-terrestrial  radiation,  atmospheric 
optical  air  mass,  and  column  length  of  precipitable  water.  If  cloud 
base  height,  zc,  and  boundary  layer  height,  h,  were  not  available, 
the  solar  transmissivity  through  clouds,  tc,  was  estimated,  using 
the  simple  algorithm, 

tc  =  (1  -  0.6FCC)   ,  (A. 1) 

where  FCC  is  fractional  cloud  coverage  of  the  sky.  If  zc  and  h  are 
available  (from  rawinsonde  data  or  aircraft  landing  reports),  the 
solar  transmissivity  for  a  non-reflective  ground  surface  can  be 
estimated  more  accurately  from  the  method  of  Liou  and  Wittman 
(1979)  .  This  method  uses  sun  angle  and  column  height  of 
precipitable  water  within  the  cloud  within  a  bivariate  polynomial 
regression  of  results  taken  from  an  accurate  multi-stream  discrete 
ordinates  model.  Currently,  we  assume  that  the  cloud  coverage  is 
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stratiform  and  confined  to  the  boundary  layer,  with  a  liquid  water 
content  of  0.78  grams/meter  of  cloud  depth.  Hence  the  only  inputs 
required  are  date,  time,  cloud  base  height,  and  boundary  layer 
depth.  The  algorithm  can  be  extended  easily  to  include  other  cloud 
types  and  water  content.  The  regression  form  is 


3    3 
i"=0  j  =  0 


(A. 2) 


where  ju0  is  solar  zenith  angle,  W  is  precipitable  water,  and  the  b;j 
are  the  coefficients  obtained  from  the  regression.  For  actual  non- 
zero surface  albedoes,  As  (default  value,  0.15) ,  we  modify  the  cloud 
solar  transmissivity ,  using  the  algorithm  of  Kamada  (1984), 


t_  = 


1  -  A. A. 


(1  -  0.12  A.)  (1  -  0.1  A,) 


(l-d0)Te  +  d<p) 


FCC 


(A. 3) 


where  Ac  is  cloud  top  albedo  (  default  value  for  stratiform  clouds 
is  0.55),  d  =  0.001068,  and  0  is  in  degrees  latitude.  The  total 
downwelling  solar  radiation,  SOLi,  is  then 


SOLi    =    (losin(0)     +    Idllf)r£ 


(A. 4) 


We  initiate  the  downwelling  thermal  radiation  computation,  using 
the  algorithm  of  Martin  and  Berdahl  (1983)  which  employs  surface 
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dewpoint  temperature,  TJp/  hour  of  the  day,  Hr,  and  pressure,  pr,  to 
estimate  the  effective  clear  sky  emissivity, 


ecs  =  0.711  +  0.0056Tdp  +  0.00073Tdp2  + 


0. 013cos(0.262Hr)  +  0.00012(pr  -  1000)   .        (A. 5) 

Tdp  is  readily  obtained  from  the  relative  humidity  and  temperature 
using  standard  formulas.  The  emissivity  is  then  modified  for  clouds 
according  to  cloud  base  height,  zc,  and  fractional  cloud  coverage, 
FCC.   Thus,  we  have 

ec  =  ecs  +  0.85FCC(1  -  ecs)  exp  ( 1 .  22X10"4  Zc)    .        (A. 6) 

Again,  boundary  layer  stratus  clouds  are  assumed  here  but  other 
cloud  types  are  readily  included.  We  compute  downwelling  thermal 
radiation  by  assuming  the  cloudy  or  clear  sky  to  be  a  grey  body 
thermal  emitter,  such  that 

IRi  =  eca64    ,  (A. 7) 

where  a  =  5.67x10*  is  the  Planck  black  body  constant.  We  obtain 
total  downwelling  radiation  at  the  earth's  surface  by  combining 
solar  and  thermal  contributions  via, 

RADi  =  (1  -  AJSOLi  +  IRI     .  (A. 8) 
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With  the  radiation  component  of  the  surface  energy  budget  computed, 
we  can  obtain  the  atmospheric  stability  and  temperature  flux  from 
the  ground  surface  to  the  air.  The  Obukhov  length  is  a  measure  of 
atmospheric  stability,  defined  as, 

L  =  -u.29/gke.    ,  (A.  9) 

where  u.  is  the  surface  layer  friction  velocity,  9  =  T(1000/pr)0285 
is  the  screen  height  (or  other  height  within  the  surface  layer) 
potential  temperature,  g  is  gravitational  acceleration,  k  is  0.4, 
the  von  Karman  constant,  and  6.  is  the  Obukhov  temperature  scale. 
The  temperature  flux  can  be  defined  as  the  statistical  correlation 
between  vertical  velocity  and  potential  temperature  perturbations, 
and  is  also  given  by 


w'e'o  =  -u.e.   .  (a.  io) 

Thus,  given  the  downwelling  radiation,  we  can  iterate  between 
estimates  of  L  and  estimates  of  surface  temperature  flux  until  both 
guantities  converged.  This  reguires  that  we  compute  both  u.  and  0.. 
u.  comes  from 


u.  =  U(z) 


Ln(z/z0)  -  ¥m  (A. 11) 


where  U(z)  is  the  mean  windspeed  at  height  z,  and  z0  is  the  surface 
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vegetative  canopy  roughness  length  (typically  "  1/7  the  mean 
vegetation  height) .  We  use  the  average  value  measured  at  the 
profile  mast,  0.07m.  We  supply  a  new  curve  fit  for  the  empirical 
surface  layer  stability  function,  ^m,  given  by 


^m  =  (1.19037  +  0.23Ln(-z/L)  )2 


(A. 12) 


which  is  computationally  more  efficient  than  previous  algorithms. 
Analogous  to  u.,  0.  is  given  by 


9.  =  <59(z) 


k 


]    Ln(z/z0)  -  *h 


(A. 13) 


From  Dyer  and  Bradley  (1932),  we  have 


^,  =  2Ln 


1  - 


1  -  14_ 

L 


(A. 14) 


To  obtain  the  ground  surface  "skin"  temperature,  0M,  we  assume  that 
the  potential  temperature  difference,  90  -  9,  between  the  roughness 
height,  z0,  and  height,  z,  is  given  by  59,  and  that  the  temperature 
difference,  90  -  9^,  across  the  laminar  layer  between  z0  and  the 
surface  is  given  by  9..  This  leads  to  the  skin  temperature 
expressed  by 


000  =  9  -  (9./k)  (Ln(z/z(1)  -  *h)  -  9. 


(A. 15) 
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Since  <56  is  not  known  initially,  0.  is  initially  set  to  zero  and  6^ 
is  initially  set  equal  to  6,  the  screen  level  potential 
temperature.  They  are  then  allowed  to  diverge  toward  equilibrium 
values  by  iteration.  The  net  radiative  budget  at  the  surface  is 
given  by 

NETRAD  =  RADi  -  e„  a  9f)04    ,  (A.  16) 

where  the  ground  surface  emissivity,  e„,  has  a  default  value  of 
0.95.  Following  the  Penman-Monteith  model  (1948,  1965), 
temperature  flux  into  the  ground  is  estimated  as,  Qg  =  0.15  NETRAD 
with  the  stipulation  that  for  stable  conditions  (L  >  0),  Q„  is  3.3 
times  larger.  The  Penman-Monteith  equation  is  used  to  estimate  the 
temperature  flux, 


-(7  (-NETRAD  +  QJ  -  w'q'0) 


w'G'q  =  -  0.84w'q'0 

(  Qcp  (Xgscc  +  7)  ) 

(A. 17) 


where  w'q'0  is  the  humidity  flux,  p  is  air  density,  cp  =  1005Jkg"1K"1 
is  the  heat  capacity  of  air,  X„  is  soil  relative  humidity,  scc  is  the 
change  rate  of  specific  humidity  with  temperature  for  saturated 
air,  and  7  =  cp/Lv  =  0.0004°K'  is  the  psychrometric  constant.  In 
turn  these  latter  parameters  are  obtained  from  standard  algorithms. 
The  Obukhov  temperature  scale  is  obtained  from  6.  =  -w'9'0/u..   In 

98 


this  scheme,  note  that  under  neutral  conditions  when  the 
temperature  flux  falls  to  zero,  9.  will  vanish  and  L  becomes 
infinite.  Thus,  to  avoid  infinities  during  iterative  numerical 
evaluation,  it  is  better  to  compute  the  inverse  Obukhov  length, 
1/L. 

From  a  rough  bivariate  analysis  of  our  results,  we  found  that  a 
useful  first  guess  is 

1/L  =  0.000674(300  -  RAD  I)  Ln ( 10  Z / Z0 ) /U2  (  z  )  .      (A. 18) 

In  order  to  cover  a  wide  range  of  radiation  values,  wind  speeds, 
temperatures,  and  roughness  lengths,  the  iteration  procedure  must 
be  quite  robust,  otherwise  convergence  is  not  obtained,  especially 
at  low  wind  speeds.  We  found  that  the  familiar  Golder  (1972) 
method  was  too  crude  to  be  useful.  Even  with  a  good  first  guess 
like  eqn.  (A. 18),  we  found  that  simple  iteration  was  unreliable, 
that  the  Newton-secant  root  finding  procedure  was  needed  for 
standard  cases,  and  that  second  order  Aitken  acceleration  was 
required  for  low  wind  speeds  and  small  roughness  lengths. 

The  Newton-secant  root  finding  algorithm  utilizes  the  form, 


-*if(Xi) 


Xi+i    x.  :   °i-t-i 


f(x.)  -  f(xM)  (A. 19) 

where  i  is  iteration  number,  x  here  is  1/L,  and  f(x)  is  the 
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difference  between  old  and  new  values  of  1/L.  The  object  of  the 
iteration  process  is  to  adjust  f(x)  to  approach  zero.  Unlike  the 
standard  Newton-Raphson  technique,  fortunately,  the  secant  method 
does  not  require  an  analytic  expression  for  the  first  derivative  of 
f (x)  which  in  our  case  is  not  obtainable.  The  Aitken  technique  is 
a  second  order  acceleration  method  found  in  most  numerical  analysis 
texts , 


X,X|  +  2     X  1+) 


Xi  +  3      Xi  _   °i-r3 


x,  +  2  -  2x,  +  1  +  x,  (A. 20) 


which  relies  on  the  Newton-secant  method  for  the  first  two 
iterations  then  computes  the  rate  of  change  of  the  convergence  from 
the  first  two  iterations  and  uses  it  to  obtain  a  refined  estimate 
at  the  third  and  subsequent  iterations. 

Once  the  value  l/L  has  converged,  we  may  use  it  to  obtain  a  final 
value  for  the  temperature  flux, 


w'0'o  =  -u.30/(gkL)    ,  (A. 21) 

from  which  we  can  obtain, 


9.  =  -w'9'o/u.    ,  (A. 22) 

and  finally  the  Deardorff  convective  scaling  velocity, 


w.  =  (gh  w'GV©)  m       •  (A. 23) 
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w.  has  meaning  only  when  1/L  is  negative,  i.e.,  the  surface  layer 
of  the  atmosphere  is  unstably  stratified,  such  that  the  air  above 
it  is  convectively  turbulent.  Turbulence  intensity  and  hence  plume 
diffusion  scales  with  w.2  under  such  conditions,  w.  depends  strongly 
on  temperature  flux  and  boundary  layer  height,  hence,  downwelling 
radiation  and  cloud  cover.  Thus,  this  series  of  equations  suggests 
strong  contrasts  in  turbulence  intensity  between  cloud  covered  and 
clear  sky  zones  in  the  coastal  boundary  layer. 

Note  however,  that  the  convective  scaling  velocity  is  based  purely 
on  surface  heating  effects  and  inversion  base  height,  without 
regard  for  turbulence  due  to  wind  shear.  This  is  a  good 
approximation  for  highly  convective  atmospheres.  However,  since 
the  edge  of  the  coastal  stratus  deck  shifts  chronically  over 
Vandenberg,  late  afternoon  conditions  with  strong  seabreeze/upslope 
winds  are  typically  only  modestly  convective  and  wind  shear  should 
not  be  ignored.  In  a  later  publication  we  will  describe  a  new  wind 
flow  and  turbulence  model  based  on  similarity  theories  extended  to 
the  outer  boundary  layer.  Therein,  we  will  discuss  our  extended 
turbulent  scaling  velocity  which  includes  the  effect  of  wind 
velocity  shear  at  both  the  surface  and  entrainment  zone. 
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Appendix  B.   Plume  Centerline  to  Cloud  Edge  Determination 

This  section  discusses  the  method  used  to  determine  cloud  edge 
positions  from  digitized  GOES-7  VAS  visible  light  images.  The 
cloud  edge  position  was  reguired  to  determine  the  distance  from  a 
given  cloud  edge  to  an  SF6  transect  plume  centerline  position  which, 
in  turn,  was  needed  to  develop  a  two  zone  dispersion  model  for 
Vandenberg  AFB. 

The  GOES-7  satellite  was  stationed  at  approximately  97  degrees  West 
during  the  time  of  the  LVDE  experiment.  Unfortunately,  the  GOES- 
West  satellite  stationed  in  the  vicinity  of  130  degrees  West  became 
non-operational  prior  to  the  LVDE  experiment  so  data  from  a  more 
obligue  viewing  angle  given  by  GOES-7  had  to  be  used.  GOES 
satellites  are  geostationary  so  their  orbital  paths  remain  fairly 
stationary  on  the  eguatorial  plane.  (GOES-7  was  located  about  two 
degrees  south  of  the  equatorial  plane  during  LVDE.) 

The  GOES-7  optical  scanning  device  we  used  was  called  a  VISSR 
Atmospheric  Sounder  (VAS) ,  an  upgrade  of  the  Visible  and  Infrared 
Spin  Scan  Radiometer  (VISSR)  flown  on  GOES  -  1,  2  and  3  (Gibson, 
1984) .  The  original  VISSR  scanned  192  microradians  per  scan  line  in 
the  visible  region  of  the  light  spectrum  -  with  8  photomultiplier 
tubes.  This  gave  24  microradians  per  pixel  in  the  N-S  direction. 
The  scan  lines  were  20%  undersampled,  implying  that  only  80%  of  the 
light  within  the  VISSR  field  of  view  (fov)  actually  impinged  on  a 
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photomultiplier  tube  in  the  array.  This  undersampling  phenomenon 
accounted  for  the  discrepancy  between  the  21  microradian  fov  quoted 
in  the  GOES  manual  and  the  geometric  (optimal)  24  microradian  fov 
given  by  dividing  the  scan  line  width  by  the  number  of 
photomultiplier  tubes.  At  any  rate,  the  VISSR  had  a  25  microradian 
fov  horizontally.  1821  scans  gave  a  full  disk  view.  The  N-S  angle 
of  the  full  disk  view  was  then: 


180      nn    _0 


1821   x  192   x  ^^   =  20.03°  (B.l) 


The  VAS  also  had  a  192  microradian  fov  per  scan  but  the  N-S  scan 
size  was  24  microradians ,  so  it  did  not  have  the  undersampling 
problem  of  the  original  VISSR.  It  also  had  a  25  microradian  E-W 
fov.  The  GOES-7  was  at  an  altitude  of  approximately  35800  km  so  the 
pixel  sizes  were: 

24  x     10"6  x     35800   =  0.859  km    (E-IV)         ,  (B.2) 

and 

21  x     10"6  x     35800   =  0.752  km    (N-S)         .         (B.3) 


The  GOES-7  VAS  data  was  acquired  from  the  SSEC  GOES  Archive  at  the 
University  of  Wisconsin.  The  data  was  sent  on  magnetic  tape  media 
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in  "GARTAPE"  format.  There  were  48  hourly  images  on  the  tape 
starting  at  1700  and  ending  at  2200  hours  daily  between  August  10 
and  August  17.  Each  of  these  images  was  examined  on  a  video  screen 
for  image  processing  quality. 

Next,  a  program  was  used  to  draw  a  40  by  40  box  around  the  domain 
of  interest  for  each  image  and  write  the  positions  and  pixel 
brightnesses  to  a  data  file  on  a  floppy  disk.  A  landmark  was  chosen 
in  each  image,  so  that  the  tracer  source  point  and  plume  centerline 
positions  could  later  be  mapped  accurately  into  this  domain. 

The  next  step  in  determining  the  distance  from  a  cloud  edge  to  a 
plume  transect  centerline  position  was  to  determine  the  pixel 
dimensions  in  the  GOES  images.  A  common  method  of  determining  pixel 
dimensions  was  to  determine  the  change  in  pixel  size  relative  to 
the  Sub-Satellite  Position  (SSP)  pixel  dimensions  using  plane 
geometry.  Several  permutations  of  this  method  were  tried  with  poor 
results  until  we  found  that  the  GARTAPE  images  as  decoded  by 
standard  processing  algorithms  was  distinctly  non-linear  and  not 
consistent  with  the  assumptions  of  plane  geometry. 

In  the  end,  the  pixel  dimensions  were  determined  locally  on  a 
visible  image  by  counting  the  number  of  pixels  from  Pt.  Conception 
to  Pt.  Sal  and  to  Purisima  Point.  These  distances  were  also 
measured  on  a  fine  scale  map  of  the  Vandenberg  region  (DMAAC, 
1976) .  The  distances  and  pixel  counts  between  these  landmarks  were 
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used  to  determine  the  x  and  y  pixel  dimensions  and  the  angle  at 
which  the  GOES  image  was  rotated  relative  to  a  local  meridian.  This 
was  achieved  by  iterating  x  and  y  coordinate  transformations  with 
an  image  rotation  angle  to  find  a  triad  of  values  giving  the  best 
fit.  The  image  rotation  was  found  as  15  degrees  (+/-  3  degrees)  and 
the  x  and  y  pixel  dimensions  1.0  and  1.3  km  (  +  /-  .1  km), 
respectively.  We  could  not  make  this  determination  more  accurately 
by  taking  more  exact  ratios  over  a  larger  domain  because  the  non- 
linear distortion  of  the  GOES  image  would  also  distort  this  ratio 
over  large  distances. 

Once  the  pixel  dimensions  in  the  Vandenberg  region  were  known  a 
program  was  written  to  determine  the  distance  from  the  cloud  edge 
to  the  plume  centerline  (pel)  position.  This  reguired  mapping  the 
source  point  and  pel  position  onto  the  GOES  image.  Source  point, 
pel  positions  and  time  intervals  for  these  two  measurements  were 
obtained  from  the  LVDE  data  set  and  encompasses  tracer  transect 
data  from  three  mobile  gas  chromatographs .  A  pixel  counting 
algorithm  was  employed  to  compute  the  distance  from  the  plume 
centerline  to  the  cloud  edge  and  also  from  the  cloud  edge  to  the 
source  point. 

For  those  pel  times  which  did  not  coincide  with  an  hourly  GOES 
image,  the  pcl-to-cloud  edge  and  cloud  edge-to-source  point 
distances  were  time  interpolated  between  preceding  and  subseguent 
on-hour  GOES  images.   Since,  in  general,  the  shape  of  the  cloud 
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front  changed  slowly  from  hour  to  hour,  we  made  no  attempt  to 
distinguish  between  distance  changes  due  to  changing  cloud  front 
shape  as  opposed  to  cloud  front  propagation. 

There  were  six  possible  sources  of  error  in  this  computation: 

1)  determining  the  pixel  location  on  the  GOES  satellite  image  which 
represented  the  location  of  a  landmark  (for  instance,  Pt. 
Conception)  .  This  error  was  random  and  was  estimated  as  +  /-  1  pixel 
vertically  and  horizontally  on  the  image. 

2)  determining  the  pixel  resolution.  The  GOES  images  were  scanned 
at  an  oblique  angle  as  well  as  distorted  by  image  processing. 

So  the  x  and  y  extent  of  each  pixel  could  not  be  determined 
accurately.  The  error  in  pixel  resolution  throughout  the  image  was 
more  dependent  on  distortion  than  on  curvature  due  to  oblique 
scanning.  We  estimated  this  error  source  at  +/-  0.2  km  throughout 
the  image  by  comparing  pixel  dimensions  using  17  different 
landmarks.  However,  it  was  closer  to  +/-  0.1  km  over  a  smaller 
domain  such  as  the  size  of  Vandenberg  AFB. 

3)  determining  the  cloud  edge  position.  The  appearance  of  clouds 
depended  on  the  angle  of  incidence  on  the  cloud  and  the  dispersion 
of  this  light  to  the  satellite  sensor.  Therefore,  the  apparent 
location  of  the  cloud  edge  varied  with  these  angles  even  though  the 
cloud  itself  did  not  change.  These  errors  could  be  partly 
corrected  in  principle,  but  were  estimated  to  be  small  compared  to 
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other  errors.  A  threshold  pixel  brightness  value  (of  84  for  the 
cloud  edge)  was  found  to  be  a  useful  value  for  all  of  the  GOES 
images.  Pixel  brightness  values  were  dimensionless  and  ranged  from 
0  to  255  for  visible  light  images. 

4)  calibrating  pixel  brightnesses.  This  problem  was  related  to  the 
error  source  in  3)  but  was  due  to  the  VAS  sensor  -  the  VAS  has  no 
onboard  calibration.  No  calibration  was  attempted  and  the  VAS  data 
was  assumed  to  be  calibrated  in  a  relative  sense  for  the  duration 
of  the  experiment. 

5)  determining  the  cloud  edge  within  a  pixel.  This,  of  course, 
could  not  be  accomplished  because  the  resolution  of  the  image  is 
limited  to  pixel  dimensions.  The  cloud  edge  was  therefore  assumed 
to  be  at  the  center  of  the  pixel  which  defined  the  edge.  The 
percentage  error  in  determining  the  pel  to  cloud  edge  distance  is 
therefore  largest  when  the  source,  pel  and  cloud  edge  are  in  close 
proximity.  For  cases  in  which  all  three  of  these  locations  are 
adjacent  or  even  collocated,  this  error  would  approach  the  pel  to 
cloud  edge  distance  itself  and  therefore  drown  out  the  resulting 
distance  determinations.  Since  the  pixels  are  basically 
rectangular,  the  maximum  value  for  this  error  is  the  distance  from 
the  center  of  a  pixel  to  a  corner  (one  for  the  pel  pixel  and  one 
for  the  source  pixel) .  For  a  pixel  of  the  dimensions  found  in  the 
GOES  images  for  this  experiment  (1.0  km  in  x  -  direction  and  1 . 3  km 
in  y  -  direction),  this  comes  to  0.8  km. 
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6)  navigating  the  true  locations  of  the  source  and  plume  centerline 
positions.  The  longitude  and  latitude  of  the  source  and  plume 
centerline  positions  were  estimated  to  be  accurate  to  +  /-  0.05 
minutes.  This  corresponds  to  approximately  +  /-  0.08  km.  This  error 
is  systematic  for  both  of  these  measurements. 
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